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Abstract 

\ These lectures start with a brief overview of sahent features of the criti- 

' cal region of hot QCD. The main emphasis is on the accurate description of 
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static plasma observables by the well-known hierarchy of reduced actions com- 
bined with 3D simulations above T ~ 2Tc. A striking pattern emerges, put in 
^ ' perspective by completing the quasi-particle picture. 

<: 

^ ■ 1 Introduction 

T-H ' This school gathers experimentalists and theorists in about equal numbers. 

. And that is not more than natural: what else is physics than a playground for 

i experiment and theory? In fact, this is the Galilean concept of physics and the 

CN ' raison d'etre of this school! 

■ So I will try to oblige and explain to you the little I came to know over the 

, last years about the plasma state of QCD. With the data from RHIC presently 

I coming out, this is an ideal subject for such a mixed audience. Also, it ties in 

• with the theoretical lectures given by Profs Strassler on flux, Shuryak on QCD 

at the transition point and Teper on lattice results. 
QhI One of the first duties of a theorist is being able to tell the experimentalist 

what to compare his data with. And the main part of this lecture will be dealing 
with exactly that: well above the critical temperature, say two to three times 
• 5^ I its estimated value of 170 MeV, a theorist can predict with good accuracy the 

^ ' behaviour of equilibrium properties of the plasma. The RHIC experimentalist 

. may be disappointed: what he sees is a far cry from a plasma in equilibrium: 

only recently monojets, increasing proton to pion ratios as function of centrality 
of the heavy ion collision have been measured and are according to some [J a 
sign of a plasma formed at about twice the critical temperature. It may be only 
at ALICE that we will attain temperatures like a GeV. At any rate it is beyond 
my competence to discuss the experimental signatures for the production of a 
plasma in equilibrium during what may be coined as a fleeting moment. In fact 
the creation of a plasma state at RHIC or ALICE is sometimes compared to a 
"little bang" . But the comparison is somewhat biased: whereas it is clear that 
the expansion of the universe was involving time scales much larger than the 
typical time scales for the QCD plasma to come to equilibrium, the time scale 
of the heavy ion colliders is down by roughly the ratio of the Planck mass and 
the pion mass. So the reconstruction of the "little bang" from the data remains 
a daunting task. 
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So as a consolation of some sort: hot QCD in equilibrium may be useful for 
the cosmologist... 

Despite recent interesting developments in high density QCD I will limit 
myself to zero density. 

In section I will discuss a very simplified version of the QCD transition. 
This to set the stage for the formal developments based on QCD. Then in 
section ^ QCD and its global symmetries and order parameters are discussed. 
Global symmetries are paramount in shaping the phase diagram of QCD. One 
of these global symmetries is explained in more detail in section (jlj). 

The change in the range of the forces [2j from the hadron phase to the 
plasma phase is the subject of section ©. The confining force in the hadron 
phase gets screened in the plasma phase. In QED it is the Coulomb force that 
gets screened. 

In QED one can introduce Dirac monopoles as external sources. One expects 
no screening at all (after all Galilei observed sun spots, witnesses of the long 
range of static magnetic fields in the sun's plasma). In QCD monopole sources 
are screened for any temperature! 

In the next section the strategy for doing perturbative calculations is set 
up. The strategy is the use of a hierarchy of effective actions. 

They are the electrostatic action, obtained from QCD by integrating out all 
hard modes of order T. 

The magnetostatic action is obtained from the electrostatic action by inte- 
grating out all fields with an electric screening mass. This leaves us with a uni- 
versal action (only depending through its coupling parameters on the original 
QCD action parameters). It is universal in its form, containing the magnetic 
field strength and covariant derivatives thereof. Six years ago Shaposhnikov 
gave a set of lectures at this school on this subject ^H], and you are vividly 
advised to read those. In the remaining part of the lecture we will concentrate 
on developments in the last four years. The main topic is the perturbation 
series as produced by the reduced actions and how well it converges when con- 
fronted with non-perturbative lattice simulations. The results seem to indicate 
that the quasi-particle picture is still a good guide, provided we are willing to 
accept not only the transverse non-static gluons but also the static transverse 
(magnetic) gluons as quasi-particles. 

In section ((T)) follow a few examples: Wilson loop, pressure, Debye screening 
mass, magnetic screening mass and 't Hooft loop. 

In section Q we come back to the gluonic quasi-particles and show how 
they determine through their fiux the behaviour of the spatial 't Hooft loop. 
In particular, simple scaling laws follow from counting rules which are backed 
by perturbative calculation. For spatial Wilson loops, analytic calculations are 
not available. We hypothesize soft transverse gluons follow analogous counting 
rules, and arrive at a scaling law for the Wilson loop which is verified by lattice 
simulations. 
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2 Qualitative picture of the transition 

To know better what we are talking about, have a look at fig.(^. 
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Figure 1: Proposed phase diagram for QCD. 2SC and CFL refer to the diquark 
condensates defined in ref . |2j . From ref . 



It is a schematic phase diagram of QCD as function of temperature and nu- 
cleon density, or more precisely the nucleon chemical potential. At the origin we 
have a groundstate where the quarks and the antiquarks combine into Cooper 
pairs. This condensate of Cooper pairs is sensitive to temperature and chem- 
ical potential changes. You see a familiar transition at zero temperature and 
chemical potential about 900 MeV: the formation of nuclear matter. For still 
higher chemical potential we get a a phenomenon called Pauli blocking. At high 
nucleon density the Pauli principle frustrates the formation of quark-antiquark 
pairs because the high density of nuclear matter renders all low lying quark 
levels occupied. So the chiral quark condensate will diminish with growing 
density. On the other hand at very high densities the gauge coupling becomes 
so small that perturbation theory is valid. It tells us the pairing of quark-quark 
pairs is preferred. Then the Cooper instability changes the groundstate into a 
state of matter where we have instead of a condensate of paired electrons with 
electric superconductivity a condensate of paired quarks with colour supercon- 
ductivity [2j. This phenomenon may take place in neutron star cores. 
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Let us now increase the temperature T. On the vertical axis at zero fi 
probably a cross-over behaviour results for realistic quark masses. Cross-over 
behaviour means a gradual change of thermodynamic quantities, like pressure 
and internal energy. Numerical simulations for realistic quark masses are not 
yet decisive on this point. Some lattice simulations of QCD [H] [7j indicate a 
critical point for non-zero nucleon density. At any rate, it is not excluded that 
the part of the phase diagram dubbed "hadronic fluid" is smoothly connected 
to the region where ALICE and RHIC will probe the phase diagram. 

The continous lines show a first order transition. First order means that 
quantities like the energy density jump at the transition. At the endpoint, 
thermodynamics tells us the transition must be second or higher order. In such 
points first derivatives of the free energy are still continuous (or even higher 
derivatives). A transition is often caused by a change in the way a global 
symmetry is realized^. We will see more about that in the next sections. 

So the diagram shows a rich variety in physics: collider experiments take 
place near the T-axis, cosmology on the T-axis, and astrophysics at low-T high 
density. That collider physics and cosmology have small density in common 
is a fortunate coincidence: one may have a direct bearing on the other and 
there is a rich litterature on this subject. From now on we concentrate on zero 
chemical potential. 



2.1 A simple model of the transition 

The simplest way to see there must be a transition is to take the bag model of 
hadrons. Increase the temperature up to energies E on the order of the pion 
mass~ lAOMeV. The Boltzmann probabilty exp— ^ for thermal excitation 
tells us a gas of relativistic pions has formed, with a Stefan-Boltzmann pressure: 

PL = 3x— (2.1) 
^ 90 ^ ^ 

There is an isospin degeneracy factor 3 in this pressure at low T. 

Similarly, coming in from temperatures T much larger than the pion mass, 
we can expect on the basis of asymptotic freedom a gas of free quarks and 
gluons. Taking the degrees of freedom into account (for a given number Nj of 
flavours) one finds: 

7 vr^ vr^ 
PH = Pgg + Pglue = 2x2x3x-Nf—T^ + 2x8x— (2.2) 

Near the critical temperature the bag pressure B of the hadrons is released 
and adds up to the pressure of the pionic gas. In other words, the individual 
hadron bags become one large bag. This is typically what percolation is about. 
Percolation of the pions in the gas means they are starting to overlap with 



"'^The Curie point in ferromagnetism is a transition where rotational symmetry is restored. Below 
the Curie point the ground state of the system is not rotationally invariant: there is a permanent 
magnetization. 
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the result that quarks and gluons do not know anymore to what hadron they 
belong. 

So, when comparing the two pressures at one finds: 

Pl{Tc) + B = ph{Tc) (2.3) 

With a bag pressure B ~ (200MeV)^ and Nf = 2 one arrives at Tc ~ 
lAQMeV. 

Here we suppose the bag constant not to vary with temperature. This means 
that the internal energy density is related to the pressure by eL(Tc) = ?>pl{Tc) 
whereas eniTc) = 2>ph{Tc)- So eq. (|2.3j) tells us that the latent heat Aec = 
^h{Tc) — ^l{T(^ = 3B, a very strong first order transition indeed! Comparing 
the jump to the mean value Cc one finds Aec/ec = 0(1). 

Such spectacular jumps would leave their marks in distributions and corre- 
lations of the hadronic decay products. 

But we mentioned already a caveat: we supposed the bag constant not to 
vary with T, and this was what made the transition first order. 

So the real question is: what does QCD say about the transition? 



3 Global symmetries, order parameters and 
the phase transition in QCD 

The QCD action has as input parameters the experimental values of A^,the 
number of colours A^, the number of flavours Nj and the masses of the 
quarks. Together with the Lagrangian: 

1 _ 
J~-QCD = -TrF^^F^y + ^ tpi{'^^D^ + mi)i)i (3.1) 
^ i=i 

these input parameters describe all of hadron physics^. 

This Lagrangian is a strongly coupled system. Its particle spectrum con- 
sists of glueballs, and quark bound states. To test this Lagrangian, numerical 
simulations with a lattice version of QCD are done. This lattice version of 
the gauge field action is in terms of SU(N) matrices U{1) living on the links 
I of the lattice. The links have length a, the cut-off in our theory. The field 
strength matrix F^y is replaced by the product of the link matrices on every 
plaquette U{P) = Ili^pU{l). This product is the exponentiated fiux through 
the plaquette: 

UiP) = Ui^pU{l) = expia^F^^ + ... (3.2) 

where we suppose the two sides of P are in the fiv direction and where the dots 
mean higher derivatives. And the action density is replaced by: 

TrF^^^^l-^ReTrUiP). (3.3) 
^Dfj, is the covariant derivative + igA^ and F^i, — d^Av — dnuA^ + ig[A^, An]. 
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The lattice coupling (5 is related to the bare coupling g hy (3 = For 
more details see the lecture notes by Prof. Teper in this volume. 

Global symmetries in QCD depend on how the fermion masses are imple- 
mented. Two extremes determine qualitatively what we know for zero nucleon 
density. All quark masses zero or all infinitely heavy. In the first case left 
handed quarks ipL and right handed quarks ipR transform under the symmetry 
group SU{Nf)LxSU{Nf)ji. The Lagrangian stays invariant, but the symmetry 
is realized in the spontaneously broken mode: the left handed quark and its 
righthanded partner do couple in the real world through a term ipL^^R, trace 
over flavour indices understood. And such a term is only invariant under a 
left handed symmetry operation in combination with the contragredient right 
handed partner. This is the diagonal subgroup. The massless Goldstone bosons 
transform as an adjoint multiplet under this group. Nature provides a non-zero 
vacuum expectation value of the left right coupling 

< i^LipR + h.c. >~ (250MeV)^. 

It transforms non-trivially under the group, so is a measure of the breaking of 
the symmetry. It is an order parameter. The Goldstone theorem assures then 
the existence of an adjoint multiplet of massless pseudoscalars. 

We have left out the two U (1) factors. One factor leaves the order parameter 
invariant and is connected to baryon number conservation. The other factor 
U{1)a transforms the condensate. But due to quantum corrections - the axial 
anomaly - it is not a symmetry of the system and the corresponding Goldstone 
boson gets a mass due to the instanton mechanism [27] . 

3.1 Universality 

The order parameter $ of many statistical systems is zero above the critical 
temperature Tc. Below Tc it is non-zero, so its behaviour is non-analytic. Should 
it jump at Tc, the transition is called first order. If it is continuous but its 
first derivative jumps, it's called second order and so on. The order of the 
transition is the same for a whole class of statistical systems and this is called 
the universality class of the transition: 

• The order of the transition is determined by the symmetry and the di- 
mensionality of the system as described by the order parameter. 

So to know the order of the transition of QCD we just take the most general 
3D action consistent with the symmetries of QCD one can write down for the 
order parameter = tpL^ii^Rj- It is the following: 

£^jj-^ = Tr{d<^'<d^) + m'^Tr^^^ + g{Tr<^''<^)'^ + hTr{'^''^f + cdet<^ + h.c (3.4) 

The critical behaviour of this action is the same as that of QCD according to 
universality. For Nj = 2 the global symmetry is that of 0(4) ~ SU{2)xSU{2) 
and is known to be 2nd order. For Nf = 3 the determinantal term drives it 
first order P^ . 
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3.2 Z(N) symmetry 

Till now the masses of the quarks were zero. Let us go to the other extreme: 
infinitely massive quarks. They leave us with only gluons as dynamical agents. 

Then there is a global symmetry: Z(N) symmetry |8j with = 3 in case 
of QCD. Z{N) stands for the subgroup of SU{N) that commutes with all 
elements of SU{N). It consists of matrices z^In, the NxN unit matrix, 
and Zk = expik^, k=0,l,...N-l. For notational convenience we'll drop the unit 
matrix henceforth. 

Where does this symmetry come from? In contrast to chiral symmetry, 
Z(N) symmetry is not a symmetry acting on quantum states. It is a symmetry 
of the free energy of the system expressed as a path integral. 

To get this point we have to understand a basic fact about the description 
of static phenomena at finite temperature T. Any observable O has a thermal 
expectation value given by the Gibbs sum: 

< O >T= Tr O exp (-H/T) /Tr exp -H/T (3.5) 

The trace is over physical states only. Physical states are by definition 
gauge invariant states, that is, invariant under gauge transformations that are 
regular in configuration space. Of course only gauge invariant observables are 
admitted. The factor 1/T in the Boltzmann factor is like an imaginary time 
span in a quantum mechanical amplitude. The transcription to a path integral 
is then straightforward [D). The trace means the path integral will be periodic 
in this imaginary ("Euclidean") time for bosons ^. 

An immediate consequence is the transcription of Feynman rules. For finite 
temperature the Feynman rules in Euclidean space undergo one single and 
simple change: instead of integration over energies, energies are now discrete 
because of the (anti)- periodicity. For bosons they equal ujn = 27rTn, for 
fermions Un = 27r(n + ^)T. In both cases n is integer. This change goes into 
the propagators, vertices and energy momentum conservation at the vertices. 

Long ago 't Hooft jSj realized that the periodicity in time does not neces- 
sarily mean you need gauge transformations to be periodic in time. A gauge 
transformation can be periodic modulo a center group element expik^ of the 
gauge group SU{N). The gluon field being an adjoint does not feel any cen- 
tergroup element. So action and measure of the path integral are insensitive to 
such a gauge transformation ^. 

So locally we have a gauge transformation. But observables that are non- 
local over the whole periodicity range will feel a change, despite the fact that as 
observables they need to be gauge invariant against everywhere regular gauge 
transformations . 

■^And antiperiodic for fermions. In this formalism boundary conditions tell an important distinc- 
tion between bosons and fermions: they tell us the distinction in statistics! 

^But quark fields are sensitive to the center group: antiperiodic boundary conditions are changed, 
and hence the statistics. 
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A prime example, where our discontinuous gauge transformation makes it- 
self felt, is the Wilson line running in the periodic time direction: 

P{Ao)ix) = ^TrVexpig dTAo(f,r) (3.6) 

where the path ordering is defined by dividing the interval [0, 1/T] into a large 
number Nr of bits of length Ar = : 

V{Ao)=\\ui U{t = 0,At)U{t = At,2At) [/(r = 1 - At, ;^) (3.7) 

in an obvious notation. We have dropped the x dependence to avoid clutter in 
the formulae. Every factor is a string bit C/(r, r + Ar) = exp igATj4o(T). Every 
string bit in this product transforms under a gauge transform as VL{t)U{t,t + 
Ar)r2'''(r + Ar) + 0{^) . So periodic gauge transforms transform the Wilson 

line like 0(t = Q)V{Aq)Q){t = 0). And so the trace is invariant. But a 
discontinuity will multiply the Wilson line with the center group phase Zk, if 

fifc(r = ^) = zl^kir = 0) (3.8) 

Note that a discontinuity other than the center group would be fixed at the 
time T = y inside the trace of the Wilson line. Only the center group is a global 
group, i.e. it does not matter where in time the discontinuity was defined to 
be. 

Although the Z(N) transformation leaves the path integral, hence the free 
energy invariant, the question whether it commutes with the Hamiltonian of 
the system makes no sense. There is no such thing as a conserved charge. 

This in contrast to canonical Z{N) symmetries that do commute with the 
Hamiltonian. They can be broken at low temperature but not at high temper- 
ature as intuition has it. In section Q we illustrate this point in QCD. 



3.3 Wilson lines, Z(N) symmetry and the deconfin- 
ing phase transition 

The thermal average of the Wilson line is related to the free energy excess AF^ 
of a state with a very heavy test quark -i/^j, averaged over all gauge transforms 
of the state and averaged over the colour indices i: 



exp-AF^/r = j Z)A^TrP(Ao)exp-5(A)/ j Z)Aexp-5(A) = {P{Aq)). 

(3.9) 

In appendix A we prove this relation. It is valid for any heavy point source 
in the fundamental representation. A source in any representation r of the 
group will just change the representation of the Wilson line into r. 
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The thermal average of the Wilson line has been simulated and it is zero 
at low temperatures, but at Tc it rises abruptly to acquire the value 1 at very 
high T. A little thought makes this clear, because of its connection to the heavy 
fundamental charge. 

The energy excess equals the energy of the fluxtube pointing from the test 
charge. As the flux cannot return to the test charge, the length of the fluxtube 
is typically the spatial size of the box. The energy equals the string tension 
times the length so is proportional to the size of the box. 

In the phase where the fluxlines are screened, this energy is finite and will 
become zero if screening is total. 

However we have swept a problem under the rug, that of the short distance 
effects on the self energy. They are still contained in the thermal average, 
ea. (|3.9|) . and contribute in terms of the lattice cut-off a: 

AF^ ~ ^. (3.10) 

For fixed temperature the lattice cut-off goes to zero exponentially fast in 
the lattice coupling. On the other hand the inverse temperature equals size 
NrU with Nr the number of lattice points in the Euclidean time direction. So 
the free energy excess due to thermal effects is to be corrected for this short 
distance effect, and to do this is in practice quite intricate IBj. 

Let us illustrate how one determines the transition, say in SU(3). You can 
ask the question: what is the distribution of expectation values of the Wilson 
line P, averaged over the space volume of our box. Mathematically one asks 
the probability E[P) of a given value P of the line to occur: 

E{P) ~ j DA6{P - P{A^ ) exp -S{A) (3.11) 

Under a gauge transform 17^ as in ea. H3.8|l the measure and the action are 
invariant. Only the line average P{Aq) picks up the factor z^, so 

E{P) = E{zlP) (3.12) 

For S'C/(3) the distribution E is shown in fig.© at the transition temper- 
ature Td- The three peaks at the center group values are equally populated, 
and the figure clearly shows invariance under multiplication with expi^. The 
central peak at P = is a sign that the system likes to be in the hadron phase 
at the same time. This suggests coexistence of the hadron and deconfined phase 
at T = Td- At higher T the central peak disappears rapidly and we are left with 
the three peaks at the center group values. This is confirmed by perturbative 
calculation of the distribution . 

So the behaviour of the Wilson line indicates that at low temperature 
the Z(3) symmetry is restored and broken at high temperature, at first sight 
counter-intuitive. It seems that at high T the Wilson line spins want to align. 
This would be understandable if the surface tension between regions where the 
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Figure 2: Thermal Wilson line histogram in the SU(3) gauge theory at the deconfining 
transition point obtained on a 24^ x 36 x 4 lattice, QCDPAX collaboration |3]. 

Wilson lines point in different Z(3) direction becomes very large at high T. 
The surface tension has dimension (mass)^. In quarkless QCD there is no scale 
so the tension must be proportional to T^. Hence alignment is energetically 
favourable at high T and the Z{3) symmetry is spontaneously broken^. 

3.4 Z(N) universality 

Let us now discuss universality in the context of the Wilson line. The Wilson 
line is an order parameter and serves therefore to define the universality class ^ . 
In analogy to the discussion of chiral symmetry, especially ea. H3.4|) . we now 
look for a 3D action which has Z(N) symmetry. The Wilson line P{Aq) is now 
written as a complex number p. In any lattice point x we have an independent 

^Note that the QCD scale A was left out of the argument. Why were we allowed to do so? At 
high T the A parameter is absorbed in the running coupling and is nowhere else present in high T 
observables (see section®). 

^For a thorough discussion see the lecture notes of Pisarski [M] . 
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"spin" and it transforms under global Z{N) as 



z^p with k = 0,1,2,..,N - 1 and z = expi— . (3.13) 

So p takes only values in the center group. 

Now it is an easy task to write down actions that are Z(N) invariant. For 
N=2 it is the famous Ising model in 3D, that models spontaneous magnetiza- 
tion: 

SN=2 = Pj2('^-p(^)p(y)) (3-14) 

n.n. 

The sum is over all links that connect two neighbouring lattice points. One 
sums the Boltzmann factor exp — S n {{p{x)}) over all configurations {p{x)} to 
get the free energy per lattice point {V is the number of lattice points): 

e^p-pVfN=2iP)= ^ e^V-SN=2{{piS)})- (3-15) 

{pm 

This gives the free energy and a transition is found at a /3c ~ 0(1). Below 
this point the order parameter (p) = 0. This is understandable: at /? = the 
relative sign of neighbouring spins does not matter in the Boltzmann factor, so 
disorder will prevail and no magnetization (p) results. 

Above this point it starts to grow to attain the value 1 or -1 at large /?. At 
large /3 the spins at the end of any link align because that lowers the action. So 
the Boltzmann probability will be higher. Whether the resulting magnetization 
is positive or negative depends on the way we prepare the system. Even an 
infinitesimal applied magnetic field h ( in the guise of a term hJ2xPi^) 
action) will decide about it. 

One can induce a region of up-magnetization next to a down-magnetization 
region by changing the interaction on the links that pierce the wall between 
the two regions. The change is from 1 — p{x)p{y) to 1 +p{x)p{y). The effect of 
that change or "twist" is for large /3 that the spins on such links will anti-align, 
because that optimizes the Boltzmann probability on such links. So this creates 
a domain wall around the twist with a surface tension a(/9) > as /9 > /?c. An 
equivalent way to create the same domain wall is to fix spins to be up at 
one end of the volume, and down at the other end. Exactly analogously, in 
gauge theory one can fix the thermal Wilson line and compute a "domain wall" 
tension. Alternatively, one can define a twist in gauge theory most naturally 
in the lattice formulation. 

The critical properties of the 3d Ising model have been well established, 
by numerical means, scries expansions etc.. The transition is known to be 
second order. So magnetization and surface tension go smoothly to zero above 
the critical point. In particular a{j3) ~ |/3 — Pd"^", with 2u = 1.26.... And 
indeed the corresponding transition in SU{2) is second order, it turns out by 
lattice simulation. In a later section we shall sec more manifestations of this 
universality for SU(2), namely for the surface tension. The surface tension 
of a wall separating two regions where the thermal Wilson line has opposite 
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signature has been simulated 07] and is shown in fig. (jllj) . Universality is well 
satisfied by the exponent. 

For SU(3) the spin action is Z(3) invariant and reads: 

5;v=3 = /?E(i-?'(^)?'*(y) + c-c-) (3-16) 

n.n. 

The first term is obviously Z(3) invariant. The reality of the action reflects 
the charge conjugation invariance of the SU(3) theory. Charge conjugation 
guarantees that the average of P and of P* is the same. And so does the 
reality of the spin action for the average of p and p* . 

Now the transition is first order for the spin system. And the SU(3) tran- 
sition is indeed first order, though weakly so |34j . With weakly is meant that 
the ratio in the jump in energy over the energy is small, in contrast to what we 
found in the bag model before. 

So universality seems to be well established for N=2 and 3. 

Not only the order of the transition but quantities like exponents are iden- 
tical according to universality. We will come back to those when discussing 
correlations. 

3.5 Universality for large number of colours 

Though QCD has three colours, GUT theories feature more, and it is therefore 
not academic to look at > 4. For N=4 and larger we find a lack of predic- 
tivity. In fact Z(N) spin models usually comprise different universality classes. 
Z(2) and Z(3) are rather the exception! 

The most general model with Z{4) invariance has two parameters instead 
of one: 

Sn=4 = Pj2(diil- p{S)p*m + c.c. + d2{l - {p{x)p*{y)f))- (3.17) 

n.n. 

The normalization of (3 is clearly a convention so we have indeed two physical 
degrees of freedom. 

The two-dimensional phase diagram of this model is well known. It has first 
and second order transitions as one varies the ratio di/d2. Setting di = the 
resulting model is Ising-like because the remaining interaction term fluctuates 
between ±1 and we saw before it was second order. Putting di = d2 = d gives 
us an action with value 2>(3d if the spins are aligned, and —(3d if otherwise. This 
is a class of models - Potts models - which has a first order transition from 
A^ = 3 on. Depending on the ratio di/d2 the order of the transition changes. 
In other words there are at least two universality classes in Z(4) spin models , 
and the question is to which one the SU(N) gauge theory belongs. 

Simulations of SU(4) [23] gauge theory show a first order transition, and 
the same is true for SU{Q) 

In conclusion: where it is defined, universality works well. For A^ > 4 we 
have to invent extra criteria to pinpoint the universality class in the spin model. 
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Figure 3: The QCD phase diagram of 3-flavour QCD with degenerate (u,d)-quark 
masses and a strange quark mass m<, jB]. 



3.6 The phase diagram of QCD 

Based on what we learnt above, fig.© indicates schematically where one can 
expect first and higher order transitions. What varies in the diagram is the 
value of the m„ = mass and the mass of the strange quark m^. The upper 
right corner contains the Z{2)) transition. As discussed before, it is known to 
be first order (2^. The deconfining transition is rather high. The lower 
left corner contains the case of the = 3 chiral limit, which is first order as 
well, according to the renormalization group analysis of eq. (|3.4() . The transition 

temperature ^ is lowest. The case N ^ = 2, upper left corner, is second 
order j26j . Its transition temperature is not as low as that of three flavours. So 
the borderline between flrst order and crossover ends in a second order point 
at niu = rrid = 0,ms = oo. Gavin et al. j25j find the critical behaviour of the 
lower left borderline is governed by an effective action with a Z(2) symmetry. 
This Z{2) symmetry is not present in the original QCD action. 

Clearly the determination of the exact location of this line vis a vis the 
physical values of the quark masses is of paramount importance. 



3.7 Chiral and Z(3) order parameters in flavoured 
QCD 

In fig.(@l) the transition region for two fiavour QCD is shown. We have only 
one symmetry, chiral symmetry. So we expect one transition at T = T^, where 
the chiral order parameter drops to zero . There are two striking observations: 
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m /T = . 06 



m /T = 0.08 



0' — L 
5.2 



Ql — ^ 
5.2 



Figure 4: Deconfinement and chiral symmetry restoration in 2-flavour QCD: Shown 
is (P) (left), which is the order parameter for deconfinement in the pure gauge hmit 
{rrii oo), and (ipip) (right), which is the order parameter for chiral symmetry break- 
ing in the chiral limit (m^ 0). Also shown are the corresponding susceptibilities as 
a function of the lattice coupling 0- 

• Despite explicit breaking of Z{3) invariance the Wilson line drops steeply 
below some T = T^. 

• The transition temperature is the same for botii order parameters: = 
Td = as peaks of the susceptibilities show. 

The first point is in seeming contradiction with the expectation that a heavy 
test quark forms easily a bound state with a dynamical light quark. However, 
you can argue that in the broken chiral phase the dynamical quarks acquire a 
mass heavy enough to recover approximate Z{3) symmetry. If so, the Wilson 
line is a sensible order parameter in the broken chiral phase. 

Confinement implies chiral symmetry breaking. After all, confinement im- 
plies a bound state of two massless quarks. But in a bound state the quarks 
must be able to flip their helicity. If so, then chiral symmetry cannot be re- 
stored before deconfinement sets in: <Ty.. In between, the Wilson line could 
be almost unity with the chiral symmetry still broken! However, Nature tells 
us that the quark condensate gets unstable above T^, and that = T^- Why 
is not understood. 

In sharp contrast, if quarks are in the adjoint representation j29jj the system 
has two exact symmetries, Z{3) and chiral symmetry. So two different transition 
temperatures are expected. Flux tubes cannot end on adjoint matter, so form 
glueballs. The region in between has no glueballs anymore, but still a fermion 
condensate and a hadron spectrum. The adjoint fermion condensate stays 
stable till T^^ ~ ST^ [IHl- 
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4 Canonical Z(N) symmetries in SU(N) gauge 
theory 

In this section we start with QCD in 3+1 dimensions. We then render one of 
the spatial directions periodic and study the effects due to varying its size. The 
results will be of later use, especially in section ((Tj). Then we will switch on the 
temperature and see what happens. 
First we fix some general notions. 

4.1 Electric and magnetic fluxes 

Below we give a quick review of electric and magnetic fluxes and their free 
energy. Let only fields with trivial N-allity couple to the gauge fields. 

4.1.1 Electric fluxes 

For example, take the y-direction periodic mod Ly and variable. The periodic 
r direction is supposed to be very long L,- = ^ — > oo, as well as x and z 
directions. 

To explain the essentials we will first put Lx = Lz = and consider the 
two-dimensional system. 

Then a time independent gauge transformation is allowed to be discon- 
tinuous mod Zk as in ea. p.Hj) but now in the periodic y-direction: 

^k{y + Ly) = zl^kiy), k = 0,l, .., iV - 1. (4.1) 
An example of such a transformation is: 

u 27r 

Ofc = exp-i-^— n (4.2) 

with Yfc = diag(/c, k,k — N, ....k — N), the diagonal k-hypercharge, an NxN 
traceless matrix with N — k entries k and k entries N — k. It has the property 
that 

(iWij;^Yk = Zk (4.3) 

is a center group element. It is a natural generalization of the familiar hyper- 
charge Yi. They span the Cartan subalgebra, which consists of all elements 
of the Lie algebra su(N) that can be diagonalized simultaneously. It is N-1 
dimensional. There is a lattice Lc of elements in this subalgebra that give upon 
exponentiation a center group element. So the elements Yfc are special points 
on this lattice. They have an important property: let q be a number between 
and 1. Then the elements Yi.{q) trace a rectilinear path in the Cartan algebra 
on which only begin and end points correspond to centergroup elements (re- 
spectively 1 and Zk)- Do this for all elements with A;=1,....N-1. Then we have an 
elementary cell of the lattice Lc- The reader can see by inspection, that this true 
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for N=3 and 4. This is the property that is useful for the Z{N) invariant Wilson 
line distribution E(P) and will be used throughout the dynamical calculations 
in subsection ((731). The stability group of Yk is SU{N - k) x SU{k) x U{1), 
the subgroup commuting with Y^. 

The gauge transformation Qk = exp iuJk is represented in Hilbert space by 
expi J dxTrDuJk-E. As operator commutes with all local gauge invariant 
operators, in particular the Hamiltonian. It does not commute with the Wilson 
line in the y-direction: 

nkP{Ay)nl = ZkP{Ay) (4.4) 

It is important that on the physical Hilbert space these operators have a unique 
effect, only depending on the discontinuity. To understand this take Vik and Q!^ 
with the same value for k. Form the product Vt^^Q.'^. In this product the singular 
behaviour drops out: both transformations belong to the same equivalence class 
through a regular gauge transformation. And a regular gauge transformation 
leaves a physical state invariant. The product of two elements from equivalence 
classes k and k' gives the equivalence class k+k' mod N. And finally $7^ is 
regular. 

As a consequence the eigenphases must be of the form expi/c^e. The 
number e is integer and conserved mod N. 

And the physical Hilbert space divides into N orthogonal subspaces TLe, 
e integer mod N, on which is diagonal with eigenvalue exp ike"^. The 
projector on such a subspace is given by 

Pe\phys >= — ^exp(-iA;e— )17fc|p/iys). (4.5) 

And since the Wilson line P{Ay) in the y-direction obeys il,kP{Ay)^l\. = 
ZkP{Ay) a state with charge e can be written as a state with e = and with 
the line P{Ay) wrapping e mod N times around the y-direction. So e is the 
promised conserved charge, and counts the number of "strings" or Wilson lines 
wrapping around the y-direction (mod N). There is a free energy for each 
of these electric flux sectors, defined by tracing only over the physical states of 
a given sector: 

exp —LrFe = Tve exp —Lt-H. (4.6) 

These free energies can be inferred from simulations on the lattice. First we 
need a formula relating the Fg to partition sums Z^, . To this end substitute 
ea. ()4.5() into ea. H4.6|) and rewrite the latter as: 

1 ^""^ 27r 
exp -LrFe = j^Y^ Z^k^ exp -ike—. (4.7) 

fc=0 

The partition functions on the r.h.s. of this equation are now Gibbs sums 
over physical states, with the operators Vik acting: 

4'^ = Trphys^k exp -H/T (4.8) 
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To understand the physical meaning of the partition functions there is an al- 
ternative definition of the operator fifc. It is only valid on the physical subspace, 
where it reads: 

An 

Qk = exp i—TrEy{yo)Yk. (4.9) 

Only on the physical subspace the two are identical. In fact they differ by 
a regular gauge transformation, as you can infer from the exercise below. 

• Show that the operator (|4.9|) has the same effect on the Wilson line P{Ay) 
as Qfc in eq. (|4.4j) . 

• Show that any regular gauge transformation Q of Qk in (|4.9j) . QQ^^"^ , has 
the same effect on P{Ay). This means that physical states stay physical 
states after applying 0,^- Hint: make use of the fact that the centergroup 
element commutes with all of the gauge group. 

So the discontinuous gauge transformation is brought about by a single 
dipole of strength at the point y = yQ. 

This generalizes easily from d=2 to d=4, adding the x and z dimensions. 
The single dipole at yo becomes a dipole sheet on the (x,z) surface at the point 
(yO)T = 0) as representing the operator (see fig.©)- We have added the 
suffix y on to distinguish it from a similar operator in z or x direction. Once 
this is done we have to admit not only operators labeled by the vector 
k = (kx, ky, kz), each component running from to N. Obviously, we also have 
fluxes e = {ex,ey,ez) and the connexion between the free energy Fg and the 
partition functions generalizes to: 

k 

exp -L.Fe- = ^ E 4'^ -^(^•^"'^) (4-10) 

k 

And so we have now a physical interpretation of the partition function eq. ()4.8() 
as the thermal average of the electric dipole sheet. It monitors the electric flux 
activity in the system as we will see in sections ((TJ and ©. 

Note that the partition functions with an electric twist are related through 

(e) 

a Fourier transform to the free energies Fg . They do not define by themselves 
a free energy as they are off-diagonal matrix elements. 

4.1.2 Magnetic fluxes 

Of course, one can define partition functions from physical states with a mag- 
netic vortex line running along a space direction, say the z direction. In con- 
tinuum language the operator creating such a vortex is: 

Vm = expiD{A)v^{x,y)E (4.11) 

with Vm{x,y) = ——^^^Y„i. When encircling the point x = y = the gauge 
transformation expivk{x,y) picks up a factor exp i ^ Yjji — Zm- This gauge 
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a) b) 



Figure 5: Electric twist a) and magnetic twist b) partition functions. The electric 
twist is a dipole sheet at fixed (r, ?/). The magnetic twist is a Dirac fiux in the z- 
direction, propagating in time, at fixed (x,y). In lattice language the plaquettes P^ ^ 
are twisted in a), P^^y in b). 



transformation remains, by definition, unchanged along the z-direction. We say 
that Vm creates a vortex or "Z(N) Dirac string". That means, a Wilsonloop W 
in the fundamental representation that encircles the vortex will pick up the 
factor: 

VmWVl = ZmW. (4.12) 

Any Wilsonloop with non-zero N-allity / will pick up a factor {zmY- But 
Z(N) neutral loops will not sense the Z(N) Dirac string, hence the name. 

Next we address the question how to propagate the string in the Euclidean 
time. As a warm up we start with a small time lapse At from r to r + At. 
Then we have for the thermal average: 

Trphys exp -HTVm exp -HAtV^ exp -H{^ - r). (4.13) 

The question is now: what happens to the Hamiltonian: 

H =^ J dx{E^ + B^). (4.14) 

The electric field strength is in the adjoint representation so so does not feel 
the Z(N) discontinuity. One might be tempted to say the same of the mag netic 
term. However, on the lattice the magnetic term is regulated as a Wilson loop 
on a plaquette. So all the plaquettes encircling the vortex pick up the Z(N) 
phase, according to eq. (|4.12|) . As the vortex runs in the z-direction the (x,y) 
plaquettes on the vortex are "twisted". 

We are interested in finite time slices. The string of twisted plaquettes in 
the Hamiltonian is then repeated in every r slice, tracing out its history. 
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The notation for the magnetic partition functions is Z)n , with the time 
shce being the full period ^: 

Zi^^ = TrphysVmew {-H/T)Vl = J DAexp -Sm{A). (4.15) 

They define directly a magnetic free energy Z^'' = exp—LT-Fjk^\ being 
diagonal matrix elements of the Hamiltonian. Note the difference with the 
electric twist partition function Z^^^ in ea. ()4.8() . 

The action 5m, is gotten from the twisted Hamiltonian. So at a given time 
slice we have along the vortex line: 

1 - ^TrU{P,,y) ^ 1 - z^{^TrU{P,,y)). (4.16) 

and repeat this for all time slices. 

The situation is shown in fig. (Ob). The vortex creates a singular dipole 
field B along the z-direction in the colour direction. 

Our definition of magnetic flux free energy has a caveat: looking at ea. (|4.15j) 
we put in a trace over all physical states. If we want no electric flux, we should 
have projected on the corresponding subspace. The reason we did not have to 
do this is the additivity of electric and magnetic flux free energy: 

Fe;rn = + F^^^ (4.17) 

Additivity is supposed to be true in the thermodynamic limit L^^y^z oo. 
Using in our definition ea. (|4.15|) the inverse of ea. (|4.1UI) we find: 

Z^^^ =Y.exp-L^Fg.^ (4.18) 

e 

For the electric twist partition function one just exchanges x and r. As the 
plaquettes Py^r deliver the electric field in the continuum limit, it is intuitively 
clear that this prescription will give the thermal average of the electric dipole 
sheet. 

4.1.3 Behaviour of flux free energies in the conflned phase 

The behaviour of the electric and magnetic free energies is quite different in 
the confined phase, where all sizes are macroscopic. 

When the size Ly of the periodic direction is macroscopic, the VEV of the 
Wilson line in the y-direction is zero. The system is confining with string 
tension a. This means that F^ — Fq = aeLy and states with e = mod 
N are energetically favoured. The energy Ef> = lim^^^^oo -^e of a state with 
e = mod N is the lowest, all others are higher by an amount ~ ct-Z^zj ^-nd 
the symmetry is "restored" because we have one unique ground state. Only 
the space 7ie=o is of importance for confining physics. It contains the glueball 
states, the localized eigenstates of the Hamiltonian. The periodicity in e mod 
N comes about because N strings decay into glueballs. 
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The magnetic free energy is decaying exponentially fast [Sj: 

- For-. aLz exp -aL^^Ly. 

This means the magnetic flux is screened. We will come back to this type of 
screening in the next section. 

4.1.4 A simple property of electric and magnetic twisted par- 
tition functions 

On the other hand the partition functions have a simple property. Suppose only 
one size L becomes small (meaning of hadronic size or smaller), whereas the 
others stay macroscopic. Consider any partition function with one single 
twist of strength k like in fig.lEt') or b)). If L corresponds to one of the directions 
orthogonal to the planes shown ( y or r in a), x or y in b)) then obeys an 
area law with the shown in the figure, and a universal coefficient pk- 

Zk ~ exp {-pkArea). (4.19) 

In the lattice formulation of the twist this universality is just a consequence of 
the Euclidian invariance under exchange of the r and x axis of fig. © a) into 
b) and vice versa. The function is computed perturbatively in subsection 
(|7.;-i|) in terms of the running coupling g(l/L). 

4.1.5 Behaviour of the partition functions in the hot phase 

Physically one would expect that electric flux free energies will show screening 
behaviour. And this is verifled easily by using the simple property of the 
partition function mentioned in the previous section at high temperature. All 
partition functions with a twist in the time direction will behave according to 
eq. ()4.19|) . These are precisely the partition functions appearing in the Z(N) 
Fourier transform that leads to the electric flux free energies in eq. (|4.2()j) : 

exp -L.f(^) = ^ exp -^{k.e''-^) (4.20) 

k 

With the partition functions on the right hand side decaying like exp — {pkArea) 
one can easily infer that the free energy differences Fg— F^ decay exponentially 
as well. 

• Deduce that i^e^o q ~ -^g^'* ~ ^^P ~Pea:^xLy 

So the electric fluxes behave radically different from the conflning phase. 
They become exponentially fast degenerate, and the electric Z(N) symmetry is 
sponateously broken. 

The magnetic flux energies behave qualitatively the same as in the low T 
phase: they are screened. 
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4.2 Breaking canonical Z(N) symmetry 

What happens as Ly becomes smaller? In fact, if Ly becomes on the order of 
the hadron size R we have a transition just as we had before with the inverse 
temperature and Ly interchanged. So we expect the Wilson line P{Ay) to 
acquire a VEV. But the physics will be that of 2+1 dimensional Yang-Mills 
plus the degrees of freedom of Ay that do not depend on y anymore, i.e. an 
adjoint Higgs. In the 2+1 dimensional space. Ay is a scalar with respect to 
rotations in the x-z plane. But there will still be a string tension a. 
So what has changed qualitatively in this phase? 

• In the x-z plane domain walls, or rather domain "lines", appear: they 
separate two regions where the spatial Wilson line, curled up in the y- 
direction, has different Z{N) values by a factor expi/c^. 

These walls have a tension pi^{M) which is perturbatively calculable with 
the 4d running coupling g{M) ( M = 1/Ly is any mass scale larger than 
the critical one ). 

The tension pk{M) is computed from the normalized twisted magnetic 
partition function = /Zq^\ At small enough Ly it behaves as 

■^fc™^ = exp —pkLxLr- This is the expression for a domainwall stretching in 
the x-direction, with energy density or tension p^^. 

What is obvious without calculation is that the tension of the wall is 

Pk = dk{l/9{M))M'' (4.21) 

from dimensional reasoning and the fact that the calculation is semi-classical. 
Its width will be 0{{g{M)M)~^). The calculation will be done in detail in 
section ((T)). 

So we find that the magnetic flux free energy in x and z direction is no 
longer screened! The domain lines are made of unscreened magnetic flux. 
This summarizes the effect of the breaking of canonical Z(N) symmetry. 

4.3 Intrinsic Z(N) symmetry in 2+1 dimensional Yang- 
Mills 

In the 3d Yang-Mills system there is also an intrinsic canonical Z{N) symmetry 
as opposed to the "extrinsic" Z(N) symmetry discussed above. It was discovered 
long ago and is due to the appearance of magnetic vortices in 2+1 Yang Mills 
theory. Hence the name "magnetic Z(N)" symmetry with as order parameter 
the vorticity. 

We give here a quick recapitulation of how this symmetry is realized and 
its relation with confinement in 2+1 dimensions |Hj ^he results are going 
to be useful in section Q. 

A vortex is created by a gauge transformation Vfc(x,z) = VLk{x^z) with a 
discontinuity expi/c^ across a line starting from the point The dis- 

continuity is not seen by the adjoint fields. Vk{x^z) has a purely local effect. 
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Only when we surround it by a Wilson loop in the fundamental representation 
it gives a phase factor to the loop: 

+ 27r 

Vk{x,z)W{C)Vk{x,z)^ = expife— H^(C) (4.22) 

if the point (x, z) is inside the loop C. 

So Vk{x^z) creates excitations that have a charge mod N with respect to 
a Wilson loop that surrounds the whole 2d system. This large Wilson loop is 
the generator of this intrinsic Z(N) symmetry. The large Wilsonloop commutes 
with the Hamiltonian. As a consequence the vorticity is conserved mod N. 

If this symmetry is realized in the spontaneously broken mode, < Vi >^ 0, 
then we have N equivalent ground states. Each of these ground states corre- 
sponds to a different orientation of the VEV and they are all mutually orthog- 
onal. 

Pick a given ground state \k > with < /i;|Vi|A: >= v Zk, < v < 1. 

Then the action of W{C) transforms it into a state W{C)\k > where inside 
the perimeter C the VEV of Vi corresponds to that in the state |A;+1 > because 
of en. (|0^ . 

• Show that Vk commutes_with |PF(C)|2 = W^QW^C). Deduce that the 
unitarized Wilson loop W{C) = \W{C)\~^W{C) has the same commuta- 
tion relation with as W{C). 

• Find < k\W{C)WiW{C)\k >= zi < k\Vi\k >, if the vortex operator acts 
inside the contour C. 

That is, the loop creates a domain "wall" [S] between the two groundstates. 

What is the typical energy density and width of this wall? In 3d Yang-Mills 
theory there is a dimensionful coupling with its square having dimension 
of mass. So a dimensional argument leads to an energy density of and the 
width l/gi- 

The VEV of the loop, < /i;|VF(C)|A; >, consists of the overlap of the state 
turned into \k + l > in the inside of the loop. So if we make the loop larger and 
larger the orthogonality of the ground states tell us the VEV goes to zero. That 
it decreases as fast as the exponent of the area follows from closer inspection 
of the overlap ^Tj. The overlap consists of a product of local overlaps j{x) 
between the vortices in a given point x inside the loop: 

< k\W{C)\k >= U^^siOli^) = exp-aS{C) (4.23) 

This argument is correct to the extent that the vortices interact only locally, 
typically over a distance defined by the dimensionful coupling g^. The tension 
(7 ~ ^3 again on dimensional grounds. 

Recall that our 2+1 dimensional system is embedded in a 3+1 dimensional 
world where the z-direction is periodic and of size M~^. The 2+1 gauge cou- 
pling expressed in terms of the 4d coupling is: 

gl = g\M)M. (4.24) 
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So for large M the gauge coupling g{M) is small. 

This means that the energy of the domain walls due to breaking of intrinsic 
Z(N) is parametrically a factor g{M)^ smaller than the tension due to the 
periodic z-direction in eq. (|4.2Hl ! And the width of the wall created by the 
Wilson loop is of the order of {g{M)'^M)^^ >> i.e. much larger than the 
extra periodic dimension. From this one concludes that, once extrinsic Z(N) is 
broken, the system is essentially 2+1 dimensional. This hierarchy of scales is 
consequential for the next subsection. 

4.4 The fate of broken Z(N) at high temperature 

Let us now heat up this 2+1 dimensional system. For the case of two colours 
the procedure is illustrated in fig.®. It shows a simulation of 3+1 dimensional 
SU(2) gauge theory with r and y directions periodic and variable, of length 
Lr = l/T and Ly = l/M respectively Pj . 

The 3+1 dimensional theory lies in region D. Simulated are the thermal 
Wilson lines Pr and the temperature Tc{M) where it becomes non-zero. The 
transition line Tc(M) starts at the M = axis at T = Tq, rises steeply and 
bends to the right. Also shown is the Wilson line Py with its critical behaviour. 
The lattice data are shown as circles and the broken lines are for our purpose 
here just fits to the data. The two broken lines are mirrored through the 
diagonal. This should be obvious: the loci of the two types of transition cannot 
distinguish between r and y: Pr{T,M) = Py{M,T). 

Start from region C on the vertical axis somewhere above M = Tq = Tc 
in fig.©. This portion of the vertical axis is the cold system in the broken 
"extrinsic" Z(N) phase. Here Py > and also the VEV of the 't Hooft vortex 
operator Vk- As we increase T along this line we first cross the broken line 
into region A where also the thermal Wilson line Pr gets a VEV. But the VEV 
< Vfc >= 0, as was argued in ref. [IHl IS]: one can show with the arguments 
of subsection ()7.3() that the correlation < Vk{x)Vk{y)^ >~ exp— m|x — y\ for 
large separation, as soon as the thermal Wilson line becomes non-zero. So this 
is the region where the intrinsic Z(N) is restored. 

But the extrinsic Z(N) is still broken because of the difference in energy 
scales. Note that the transition of the essentially 2+1 dimensional system is at 
a /izf/Zier temperature than the 3+1 dimensional system. This is intuitively rea- 
sonable. After all, in 1+1 dimensions the transition is at infinite temperature. 

Finally increasing T even more, the restoration of the extrinsic Z(N) takes 
place at the second crossing of the horizontal line into phase B. Then we are in 
the high temperature phase B of 3+1 dimensional QCD. 

The symmetry in the figure between phase B and C is deceptive from a 
physics point of view. In phase C there are one dimensional "domain" walls 
tracing out a two dimensional sheet in time r. They do separate regions where 
the Wilson lines Py have different center group values and an observer in the 
(x,z) world can observe those walls with their high energy density. Phase B is 
the hot QCD phase where there are regions with different center group values 



23 



100 k 



lOkE 



1 ki 



M/T„ 



100 



10 



D 



O' 



B 



uL 







A: IP l>0, IP l>0 

X y 

B: IP l>0, IP 1=0 

T y 

C: IP 1=0, IP l>0 

X y 

D: IP 1=0, IP 1=0 



uL 



10 



100 
T/T 



1 k 







10 k 100 k 



Figure 6: The 4d phase diagram in units of Tq = Tc, as shown by Montecarlo data, 
ref. jini- The data show the points where the Wilson hues and Py undergo the 
transition. Note the symmetry t ^ y. 



for P-r, separated by two dimensional sheets as discussed at the end of subsec- 
tion ()3.3|1 . These sheets cannot be interpreted as domain walls because they do 
not extend in time as solitons. In section Q their role will be seen to be that 
of detectors of electric flux. The tension is indeed symmetric with respect 
to the diagonal. 

The reader may ask the justified question: how do the currently popular 
extra periodic space dimensions of size 1 /TeV to our four dimensional world fit 
into this picture? The answer is simple and expected from the one dimensional 
domain walls found above: for 5D SU(N) gauge theory, phase C now contains 
two-dimensional domain wallsl The energy density of these walls is typically 
(TeV)'^. The form of the phase diagram is qualitatively the same. Hence the 
extrinsic Z(N) symmetry gets restored at high T and the walls put constraints 
on cosmological models [SU] [1^ . 



5 Forces and screening in the plasma 

Till now we did not mention what happens to the forces in the QCD plasma. 

In this section we will first discuss Debye screening on a perturbative basis. 
IN QCD this turns out to be insufficient. We need a non-perturbative definition. 
Then, in the next subsection an operator formalism is presented, useful for the 
systematics of the lattice calculations. 
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Finally magnetic screening is defined. This is a new and important aspect 
of thermal QCD. 



5.1 Electric screening 

In a QED plasma one would like to know what happens to the electric field 
due to a heavy charge Q and to the Coulomb force between two static charges 
of opposite signature at distance r in the z-direction. Let us look at scalar 
QED, as it shows in leading order in the coupling some features in common 
with QCD. 

The Coulomb force is transmitted by the the Aq potential. The propagator 
of Aq is renormalized by the one loop scalar and seagull diagram and gives the 
self energy, as calculated by the Feynman rules discussed below eq. (|3.5|) : 

nMrt-^TE/^fP^^^Mpi^.^if) ,5.1) 
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The self energy is transverse, p^Il^i,{p) = 0. It has two independent tensors, 
that we choose to be IIoo and H^/^. For T = they are proportional. 

For the Coulomb force we are interested in the static limit po = 0. We 
resum all the self energy bubbles to get the propagator and the static part of 
the < AqAq > propagator becomes : 

+ Uoo{po = 0,p) ^^''^^ 

To find the pole to lowest order in the coupling, we let p — > and find 
noo(0, 0) = We use dimensional regularization so that the 1^ = 

contribution to IIoo is ~ J dl^ = 0. Only hard modes proportional to T inside 
the loop contribute to the pole mass! 

In configuration space this leads to Coulomb screening: 

1 1 , e2T2 ^ ^ 

> — exp —rriDr with mu = . (5.3) 

r r 3 

In scalar QED the self energy is gauge independent to all orders in perturbation 
theory. And so you can start to evaluate the corrections to the Debye mass by 
computing the corrections to the pole location. Indeed, one can reformulate 
perturbation theory by adding the screening mass term to the action of scalar 
QED, and use the Feynman rules with the modified static propagator above: 

Sqed = {Sqed + ml [ Al} - ml [ a^ (5.4) 

Jx J X 

To avoid double counting you have to subtract the screening term as well, 
and use it as an insertion whenever you have a self energy sub diagram. This 
procedure leads to a well-defined perturbation series. However the powerlike 
infrared divergencies are now cut-off at ~ eT so we can expect terms in 
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the series like e'^/rriD ~ e^, i.e. odd powers in the couphng. But otherwise the 
series can be computed to arbitrary order by taking the Debye screening into 
account. 

The non-abeUan case to lowest order is qualitatively the same. The Debye 
mass changes only by replacing — > g'^{N + Nf/2). The salient differences 
are: 

• In the non-Abelian case the self energy is not gauge independent. For the 
pole location one can argue that it is gauge independent. 

• More important, already in next to leading order an infinity of diagrams 
in the static sector contributes ^Hl !30j . 

This means that an observable is needed to define the screening mass inde- 
pendently of perturbation theory. 

Indeed, there exists a natural definition in terms of the correlation of two 
static charges in the fundamental representation. From the results of the previ- 
ous sections and appendix A we can write this as the correlator of two Wilson 
lines: 

< P(r)P(0)t >= exp -^^^ (5.5) 

This can be simulated on the lattice by non-perturbative means. To this 
end one takes a lattice periodic in all directions, and the Wilson lines separated 
over a distance r in the z-direction. 

In the confining phase, below Tc, this correlator falls off at long distances 
as exp— ^^^r, due to the string tension cr(T) = VT{r)/r. Above Tc the string 
tension gets screened by the Debye mass and the potential becomes: 

FE{r) = Feo - — exp -mor. (5.6) 

r 

The parameters in this free energy are only depending on T. 

5.2 An operator formalism as a bookkeeping device 

The path integral of the spatial correlator can be read in an alternative way. 
Consider the fictitious Yang-Mills Hamiltonian in the (x, y, r) space, with its 
physical Hilbert space. This Hilbert space will contain eigenstates of the Hamil- 
tonian which are different from the one in which we live. One space direction, 
r, is finite and periodic with period 1/T. That means that the rotation group 
in these three dimensions is reduced to SO{2) times the discrete rotation group 
admitted by the periodic finite r direction. The path integral reads in terms of 
this Hamiltonian H and the said physical Hilbert space: 

< P(r)P(0)t >= {Trphysew -HL,P exp -HrPU^p -HL,)/{Trphysexp -HL,). 

(5.7) 
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The Wilson lines P are now expressed in terms of the canonical operators Aq, 
Lz = {Lz — r)/2, and Lz >> r. In the limit of ^ cxd the correlation becomes 
a some of exponential decays: 



exp —ninr. 



(5.8) 



n 



The mass gap m„ is the value of the energy compared to the groundstate energy, 
at zero momentum Px = Py = Pt = ^■ 

We want an efficient bookkeeping system for the states excited by the Wilson 
line (and eventual other observables). 

To this end we define in our fictitious Hilbert space a parity transform P, 
under which only the y-direction changes sign (and hence only Ey and Ay). 
Remember the rotation group is reduced to SO{2), so simultaneous flipping of 
X and 2/ is a rotation. 

Charge conjugation is as usual: — > —A^^. 

Still another quantum number is i?-parity:it changes r into y — r, and Aq 
into —Aq^. 

So the symmetry group is 5*0(2) x Z2{R) x Z2{P) x Z2{C), hence states 
labelled by J^*^- The 5*0(2) group is generated by xdy — ydx- So PJP = — J. 
Look at any eigenstate |j) of J with j / 0. Then the state P\j) has J = — j, 
so is orthogonal to Prom these orthogonal states we can form the parity 
doublet, degenerate in energy: 



Por spin zero states this argument fails. And indeed, lattice simulations 
reveal j46j differences up to a factor two in spin zero parity conjugates. 

Clearly our Wilson line operator excites spin zero and positive parity states. 
Also, under C and R: 



An important caveat is due to Arnold and Yaffe [41j : the potential consists of 
two channels of exponential decays! One is governed by the correlation of ImP 
and the other by that of ReP. The first is odd under charge conjugation [Aq — > 
—Aq), the second even. So they do not mix. The Debye mass corresponds to 
the odd channel, as it should, according to our definition in terms of self-energy. 

There is no difference between the two channels if < PP >=< P^P^ > = 
0. And in the phase where Z(N) is unbroken it is not hard to see that both are 
exponentially small with respect to the correlation < P^P >. The exponent is 
controlled by the string tension a. In that case the lowest energy state with 

^The combination of time reversal T {Aq AqjAi —Aj) and charge conjugation C in the 
usual Euclidean version of the theory gives an operation R = TC that only flips the sign of ^o- 
Time reversal has no effect on the Wilsonline, bwcause it inverts the time ordering and at the same 
time transposes Aq as a colour matrix. So the Wilson line is transposed as a matrix , but its trace 
stays the same. So R has the same effect on the Wilson line as C alone. 



|j;±) = (l±P)|j). 



(5.9) 



P{Ao) ^ PiAoV. 



(5.10) 
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energy Eq{T) has e = (ea;,ej^,e^) = 0. All states with eV / o are exponentially 
suppressed by factors exp —aL^, and the reader can convince himself by going 
through the exercise below that both correlators give the same area law in 
the hadron phase. This is intuitively expected from two like charges being 
unscreened: in a periodic volume their expectation value should be zero. 

• Show < PP >= ^rp/,y^e-(^-^o)^-Pe-(^-^o)'■Pe-(^-■^o)^^ up to expo- 
nentially small terms, if Z{N) is unbroken. = {Lz — r)/2. 

• From the above, show that < PP > is a superposition of amplitudes 
involving the flux states le,- > by using the projection operator ea. ()4.5() . 

• Find the suppression factors in front of < e,- + 2\Pe~^^~^°^^P\er > for 
all Cr and their absence in front of < 0|P"''e~*^'^~'^°)^P|0 >. 

In the hot phase there is no reason the two channels are the same. Two like 
charges with screening can coexist in a periodic volume. 

We can abstract the following conclusion from the above. The fictitious 
Hamiltonian has for T below Tc (don't forget that this means for H that the 
spatial dimension in the r-direction gets large enough, wheras the temporal 
direction in the z-direction stays infinitely long) a Z{N) symmetry which is 
the canonical symmetry discussed in section As long as the symmetry is 
unbroken (for T < Tc) this Hamiltonian has winding states, labeled by the 
conserved quantum number e,-. The winding states have a tension a. In the 
limit that T = we have the usual 3d Hamiltonian with the glueball mass 
spectrum and with ground state | Cr = 0). All other winding states have 
infinite energy in this limit. For finite T the winding states have finite energy, 
and are the groundstates of identical towers of states. 

On the other hand, when T > Tc, we have Z(N) realized in the broken mode. 
The N winding states |e,-) have become the degenerate groundstates with, 
again, each an identical mass spectrum. The Debye mass is one of those mass 
levels, as we will see in the next section. 

5.3 Screening of heavy magnetic charges 

Not only the force law between heavy electric charges like the heavy quark, but 
also the force between heavy magnetic charges tells us about the medium. The 
original idea of 't Hooft and Mandelstam was that of a dual superconductor, 
with the electric Cooper pairs replaced by some form of magnetic condensate. 
Especially the lattice community has been fascinated through the last 25 years 
with this idea because it defies perturbative access. 

In section (|4.1.2j) we constructed an operator Vk creating a magnetic flux of 
strength expifc^, ea. ()4.11() . 

To get the monopole anti-monopole pair at points (0, r) we have the vortex 
end at and r on the positive z-axis. The vortex is given by a gauge transfor- 
mation Vk{x) which is discontinuous modulo a center group element expik^ 
when going around the vortex. The vortex is like the Dirac string in QED. It is 
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Figure 7: Monopole antimonopole pair induced by twisting the plaquettes pierced by 
the Dirac string. 



unobservable by scattering with particles in the adjoint representation, as long 
as it has center group strength. 

The Gibbs trace can be worked into a path integral along the same lines as 
in section (|4.1.2|) . and on the lattice it takes the form |45j : 



exp 



-FM{r)/T= (^j D^exp-%)(.4)^/ j DAeicp-S{A). (5.11) 

The action S'^^,) is the usual action, except for those plaquettes pierced by 
the Dirac string. Those plaquettes are multiplied by a factor expzfc^, as in 
fig.©. 

• Show that any deformation of the string can be obtained by a change of 
integration (link) variable through a centergroup element. 

The reader will recognize from section Q the magnetic vortex but now with 
endpoints, where the monopole pair resides. Varying the endpoints permits one 
to find the potential for all separations. 

Screening is expected in both confined and deconfined phases: 



FM{r) = Fmo - CM 



exp —iriMr 



(5.12) 



All parameters are function of T. In the cold phase the screening is a conse- 
quence of the electric flux confinement. This is natural because the ground 
state contains a condensate of "magnetic Cooper pairs" , according to the dual 
superconductor analogy. It is a screening mechanism whose details are not un- 
derstood^. We dropped for notational reason the dependence on the strength 
k of the monopole. It is important to note that this strength comes in periodic 
modulo N! So the screening length is a periodic function of k and, because of 
charge conjugation, of N — k. It would be interesting to see what this depen- 
dence is. 

In the hot phase there are indications from spatial Wilson loop simulations 
that there is additional thermal screening from magnetic quasi-particles, as 
discussed in section (jHJ. 

Analogous to the Wilson line correlator we consider the Hamiltonian H in 
the fictions system of (x, y, r) variables. We search the operator Vk acting on 

^Nevertheless a quantitative understanding of the energy of a magnetic flux exists |3 as mentioned 
in subsection H4.1|l . 
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the Hilbert space of physical states of this Hamiltonian, that reproduces the 
path integral ea. (|5.11|) ^ . So should create a vortex in the (x,y) plane 
at every time slice r and the Hamiltonian H should propagate every one of 
these vortices in the z-direction over a distance r. So Vk is the 't Hooft vortex 
operator discussed in section @: 



with Vk{x,y) = arctan{l^Yk. 

Both under parity and charge conjugation the vortex Vk transforms into 
Vl Its spin J equals 0, despite the appearance of the rotated singularity 
line. On physical states the location of the singularity does not matter. Hence 
the operator ImVk excites spin zero states with P = C = —1. The mag- 
netic screening mass should correspond to the self-energy of a magnetic gluon, 
just like the correlator of the thermal Wilson line had to correspond to the 
self energy of a temporal gluon, So we choose the negative charge conjugation 
component ImVk- The magnetic screening mass distinguishes itself from the 
electric screening mass by the opposite parity. This will prove important! 

For SU{2) gauge theory IniVi = and we have to take ReVi exciting 
J^'-^ = 0"*"+ states (All of the physical Hilbert space is C = +1.). 

• Convince yourself of the quantum number assignments (parity in two di- 
mensions flips only one direction). Prove that R = TC parity is +1. 

Perturbation theory is not reliable and we need lattice simulations |3H1 |47j . 



6 A quantitative method: perturbation the- 
ory and dimensional reduction 



Lattice simulations are our only tool today for tackling the critical region of 
QCD in a quantitative fashion, as far as the problematic fermions with small 
(including realistic) masses are avoided. But the region above a few times 
the critical temperature can be accessed by the method of dimensional reduc- 
tion without that problem. As we will see, fermions come in through the 
parameters of the reduced theory. 

The method of dimensional reduction permits one to do perturbation the- 
ory, not only at very high temperatures but down to T ~ 2Tc. To obtain all 
coefficients of the perturbation series, one has to do dimensionally reduced lat- 
tice simulations, i.e. simulations in three dimensions. This is due to the three 
dimensional magnetic sector of the theory being a confining theory. 

In fact the idea is similar to that in Kaluza-Klein theories: at high tem- 
perature the periodic dimension is very small with respect to the typical mass 

^We use the same notation as for the vortex operator in {x, y, z) space as there is no risk 
confusion. 




(5.13) 
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scale of 3D Yang-Mills, and the Fourier modes in the periodic direction are 
proportional to lirTn with n integer: 



Po 27rnT, n integer for bosons. 



(6.1) 



For the fermions the Fourier modes are anti-periodic so the n are half-integer. 
All non-zero modes are called hard modes. 

So one works in three dimensional space with a 3D action. The field variables 
are the constant modes A^{x). The parameters in that action take the effects 
of the temperature into account. Some of the constant modes, the electric ones, 
get a mass due to Debye screening. The magnetic modes stay massless, at least 
in perturbation theory. So if one is not interested in distances on the order of 
the Debye screening but in much longer distances, integration of the electric 
screening modes is mandatory. We are then left with a 3D theory with only 
magnetic modes. They interact with a dimensionful coupling gj^j and describe 
a theory which is accurate on mass scales equal to or smaller than gj^j. 

6.1 Integrating out the hard modes 

To be precise we want to integrate out all degrees of freedom in the original 
QCD action that relate to momenta and Fourier modes of order T. So we need 
to fix a cut-off A^; somewhere in between the scale T and gT. In low order we 
can do without. This is because we are interested in amplitudes with external 
legs with n = and p ~ gT. To one loop order all modes with n 7^ in the 
loop introduce a scale of order T in the loop integration over p. The mode 
with n = has momenta of order gT injected from the external legs, so the 
momentum integration will involve only gT 

The form of our effective action Se is dictated by all symmetries, global and 
local of the original QCD action and which are respected by the integration 
process. That implies all the symmetries we knew already, except that the 
electric term in the static action will have no OqA term. So Aq appears as 
an adjoint Higgs term in our 3D gauge theory. The electrostatic QCD action 
density reads: 



Because of R- conjugation invariance (Aq —Aq) the electrostatic action 
must be even in Aq. 

So far for the form of the action. The parameters in the action are all 
expressed in even powers of the QCD coupling. That is because only hard 

^°To two loop order there can be internal propagators with n = and momenta on the order of T. 
To compute to this accuracy one either introduces the cut-off we just discussed or one exploits 
E{P) as generating functional dl for the electrostatic action. 



= Tr{D{A)AQf + mlTrAl + XE{Tr{Al)f + 
+ \E{Tr{AQf - \{TrAlf) + TrPf^ + 5Ce. 



(6.2) 
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modes are present in the integrals. Odd powers will appear as soon as we 
admit modes of order gT. 

The parameters in the Aq sector are needed up to two loop accuracy ^S] 
and we give the result for A^c = 3 and Nf. 

ml = gif^fil + ^)T\l + A5^(^)) and = ^r(l - ^)(1 + 6g\f^)). 

(6.3) 

For A^; see reference The coefficients A and d depend logarithmically 
on the scale // and for their explicit form see refs. jl^l - 

The gauge coupling starts to run and in the MS scheme one finds [T5] : 

gl = g\fi)T{l + log - — i- log -^)}. (6.4) 

The parameter = 47rrexp(— 7 — ^) = 6.742..T follows from the one 
loop renormalization of the F^j term through the effects of scale T ■ Eulers 
constant 7 equals 0.577214... If one subtracts at this scale the renormalization 
effects appear only to two loop order and the coupling is then function of 

.2 AT 0/1^2 



g%N _ 247r 



T lllog(^)- ^^-^^ 

MS 

Quenched QCD lattice simulations give us the critical temperature in terms 
of the QCD scale 13^1: 

1.15 ± 0.05, = 3 and — = 1.23 ± 0.11, N = 2 (6.6) 



From the three dimensionful quantities in this Lagrangian we can form two 
dimensionless quantities: 

A 171% 

X = and y = —x- (6.7) 

Te 9e 

We have made our promise true that the scale A-p^ only goes in through 
the running of the coupling. 

The dimensionless couplings x and y contain both g'^{iJi). Eliminating the 
latter gives a very simple relationship between the former: 

xy\iD = ^[l+'^x + 0{x^))ioi N = 2 (6.8) 

and 

xy|4D = ^(l + ^x + 0(x2))forAf = 3. (6.9) 

This is the trajectory in the x-y plane of the 4D physics, to order O(x^). 
We put Nf = 0. Remarkable is that it does not depend on the subtraction 
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scale /i! The subtraction scale survives of course in the variable x but not in 
the relation between x and y. If this trend continues in higher orders, the series 
in X is probably well convergent. 

Of course the physics of this effective action is specific to the quark-gluon 
plasma. First the coupling should be sufficiently small, and the presence of 
the mass term indicates that the electric flux is screened. And indeed m^; is 
to lowest order identical to the Debye mass niD since both equal the one loop 
static self energy at zero momentum. The difference comes in the corrections. 
Whereas the corrections to m^; are 0{g^) due to the absence of soft modes, 
those to ruD are 0{g) due to presence of soft modes. The Debye mass m^i is a 
physical quantity, is a parameter in the electrostatic action. 

The mixed sector is known to one loop Jl] up to six external legs. We 
lumped it into the term 5C e ■ The reason for doing so is a question of accuracy. 
Already the superrenormalizable terms retained in ea. (|6.2() do insure that the 
error we make in calculating some observable O with our electrostatic action is 
0{g^). This error to be numerically small constitutes one of several constraints 
on the value of the coupling. It warrants the calculation of the mass and four 
point coupling to two loop order above. 

Let's see how this accuracy comes about. The argument is dimensional and 
based on the invariances of the reduced theory. These are the discrete spatial 
symmetries, 3D rotational and gauge invariance. Including two extra spatial 
potentials on any of the terms va. Ce you get six independent terms JTl'l ( from 
F^ {DFf.AlF"^, AqFDF). A typical term reads: 

5Ce ~ ^{DFf. (6.10) 

The square of the coupling appears because of the interaction of the stationary 
modes with the heavy modes. The scale T is there for dimensional reasons. 
The question is what this vertex is going to contribute. Irrespective of the 
observable in which it appears, we can say that the covariant derivative D 
concerns momenta in the effective theory of 0{gT). That gives an g'^ factor 
in front of the F"^ factor already present in the original Lagrangian Ce, and 
provides the order of the relative error. 

The estimate is generic. It can be higher order for specific observables. It 
motivates the two loop accuracy for tue and A^; above. 

Another constraint is the following. The cut-off for our theory is 27rT. The 
Debye mass (~ gT) , a typical scale of our electrostatic theory, should then be 
smaller than this cut-off. This means < 27r, or < vr. 

In terms of temperature scales: if -r ^— ~ 2 then our coupling through 
ea. (|6.5j) with = 3 equals g ~ 1.7, consistent with the cut-off limit. 

6.2 Integrating out the electric screening scales 

On the other hand there is the scale g^ = g^T. If this scale is much smaller 
than the Debye mass, we can integrate out the scale gT by integrating out 
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the Higgs field Aq from Se- This necessitates the introduction of yet another 
cut-off Am separating the scales gT from g^T. That will lead to a new action 
Sm with only the magnetic fields present. This magnetostatic action density 
reads: 

Cm = TrFf^ + 6Cm (6.11) 

with a magnetostatic gauge coupling g^^ . 

This coupling is related to the electrostatic coupling gE through the renor- 
malization of the magnetic gluon field strength: 

4^(1+51^)4 

To one loop order the Aq field is the only field contributing. To compute the Z 
factor we have to compute a diagram like in fig. (|l()d ) with the wavy external 
legs the background magnetic potential with momentum of 0{g'^T). There is 
also a tadpole-like diagram contributing. 

Simple power counting gives a linear infrared divergence for the transverse 
result Z. Since the infrared in Se is cut off by niE, we expect parametrically 
Z ^ g'^/ruE g- So odd powers of g are to be expected. For gM we get for 
three colours: 

»« = ^^(l-^) (6.12) 

for all reasonable couplings g < 1 a small effect. 

Using the magnetostatic action at scales g^T or smaller will induce an error 
0{g^) with respect to the results one would have got with the electrostatic 
action. Like in the previous subsection a generic estimate tells us for a typical 
term from the correction term in eo. 1)6. 11(1 : 

5Cm ~ -^ADFf (6.13) 
ruE 

The coupling g^ describes the interaction between electric and magnetic modes, 
and niE the scale of the integrated degree of freedom ^o- Now D ~ g'^T and 
the relative error is 0{g^). 

This magnetic action has no dimensionless couplings like the electrostatic 
one. At scales gf^^ it's obvious we cannot form a small dimensionless number 
with the coupling gj^j. So the coupling in this theory is strong. A formal pertur- 
bation expansion of say the free energy gives from four loop order {0{g^)) on 
powerlike infrared divergencies as naive power counting shows. Regulate with 
a mass m. Now, any free energy diagram with L loops has a power (f?!/)^'^"^^ 
in front of the integral. So the integral must give a result m^^~^') to get the 
correct dimension for the the free energy. For L = 4 one expects a logarithm in 
the ratio cut-off over mass, calculated recently j20j . For L > 4 we have linear 
or higher divergencies. For a superrenormalizable theory all logs containing the 
cut-off are contained in L = 4. 

When one regulates these divergencies with a mass of 0{g\j) higher loop 
diagrams are all of order g^j^ modulo logarithms. This is Linde's argument [T9] . 
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What this means is that the coefficient of the sixth order free energy is not 
perturbatively calculable. We need non-perturbative input like the lattice. 

If one needs higher order effects then the term 6Cm has to be expanded in 
the magnetostatic partitionfunction Zm'- 



j DAexp ( - Sm{A) - J cLxSCm) (6.14) 
= j Dlexp -Sm{A){1 - j dx5CM + ) (6.15) 

This gives an expansion for —y log Zm = ag^T^{l + bg + ....) where a, b 
and higher order coefficients have to be computed on the lattice or any other 
non-perturbative method. 



7 Dimensional reduction at work 

We will treat a few examples of relevant observables in order of mounting 
complication. 

We will start with the spatial Wilson loop. Then the magnetic screening 
length niM- Then its electric analogue, the Debye mass and finally the spatial 't 
Hooft loop and the pressure. We want to calculate the first terms in the series 
up to and including the term where the magnetostatic action enters for the 
first time. We just saw that this coefficient has to be computed from 3D lattice 
simulations. It turns out to dominate for any reasonable temperature, say from 
a few times Tc till lO^Tc. Then one compares to 4D lattice data to determine 
the remnant of the series. This remnant turns out to be small, typically on 
the order of 30% up to T ~ 2Tc at least for Wilson loop and Debye mass. For 
pressure, magnetic screening length and 't Hooft loop, this program is being 
pursued. 



7.1 Spatial Wilson loop and magnetic screening mass: 
a window on the magnetic sector 

The spatial Wilson loop is given in terms of a spatial loop L, and a represen- 
tation r of SU{N). The vector potential in this representation, Ar, appears in 
the loop as: 

Wr{L) = TrVexpi j gAr.dl. (7.1) 

As for the Wilson line, eq. (|3.6|) . the exponential is path ordered and hence 
invariant against regular gauge transformations . This spatial loop should mea- 
sure the magnetic flux of any fixed gauge field configuration, as suggested by 
its abelian analogue and Stokes theorem. There is a useful version of Stokes 
theorem for the non- Abelian case |l2] |43j . 
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The thermal average of the spatial Wilson loop shows area behaviour with 
a surface tension ar{T). 



{Wr{L))T = exp -ar(T)A{L) + . . . , 



(7.2) 



The dots indicate perimeter terms. As far as the tension is concerned, it 
is very plausible that it only depends on the number of quark minus the 
number of anti-quark representations constituting the representation r. This 
number is called the N-allity k of the loop. Also the tension is periodic in k 
modulo N. For the tension of the 't Hooft loop these properties are verified 
easily from its definition. 

A useful corollary: a loop with N-allity k ( a k-loop) and a — A;-loop have 
by charge conjugation the same tension. So because of the periodicity also the 
N — k loop has the same tension. 

For = 2, 3 only one tension results, because of charge conjugation and 
periodicity. 
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Figure 8: The temperature over the square root of the spatial string tension versus 
T/Tc for SU{?)). The dashed line shows a fit according to a two loop scaling formula 
for the coupling, see text below eq. ()7.5p . From ref.j 



We now establish that the loop is a perfect window on the magnetic sector. 
Its thermal average in path integral language is: 



{Wk{L)) 



DAoDAWk{L)ex:p-S{A)/ / DAQDAex^]i-S{A). (7.3) 



^See the notes of Prof. Teper in this volume. 
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Integrate out all hard modes. They will not contribute to the tension 
of the loop, because the tension is due to correlations of the potential A over 
macroscopic distances. That will leave us with S replaced by Se on the r.h.s. 
of the average. Since the spatial loop contains only spatial potentials A we 
can integrate over to obtain Sm from Se- We arrive for the tension at the 
average: 

exp -ak{T)A{L) = J DAWk{L) exp -Sm{A)/ J DA) exp -Sm{A). (7.4) 

The hard and electrostatic free energies fh and /e drop out in the ratio. 

The only dimensionful scale in the magnetostatic action is g^^. So the 
tension, having dimension (mass)^, can be written as: 

Tk{T) = Ckgij{l + 0{g')). (7.5) 

So the dominant contribution to the tension is entirely from the magneto- 
static sector. In figure ( |H1) you see a fit of the tension data to this parametric 
expression for SU(3). The authors took for the magnetic coupling = g'^, 
so neglected renormalization effects of the scale gT, which are a few percent at 
T = 2Tc, see ea. (|6.1'2p . On the other hand they included two loop renormal- 
ization effects. Dropping those effects, and taking into account the uncertainty 
in the relation between A.j^ and Tc there is still consistency between data and 
the one loop formula eq. (|6.5j) . 

Notably the value of the tension at the critical temperature is within errors 
equal to the tension at zero temperature. So the tension of the spatial Wilson 
loop does not change within errors in the hadron phase. 

The conclusion is quite clear: down to temperatures a few times Tc, the loop 
behaviour is determined by leading order magnetic sector effects! These effects 
are embodied in the dimensionless number Ck=i- The rest of the T-dependence 
is through the hard- mode-running of the coupling, eq. (|6.5|) . The number Ck=i 
is within errors equal to the purely 3D simulation of the loop. 

The spatial Wilson loop measures in a sense to be specified later the mag- 
netic flux in the system. The tension is flat from T = to T = Tc, according 
to the data. In all of the conflned phase the magnetic activity does not change. 

Above Tc it starts to grow like gf^T'^. Apparently beyond the transition the 
activity goes up, and comes, as the data tell us, entirely from the magnetostatic 
sector. 

This window on the magnetic sector spurs an obvious question ^7] : what is 
the dependence of the coefficient on k and A^? The lattice data by Teper and 
Lucini j23 Siie consistent within a percent to the simplest possible picture for 
the 3D magnetic sector: that of a gas of almost free quasi-particles, in some 
sense "static transverse gluons". We'll come back to this in the last section. 

7.1.1 The magnetic screening mass 

The magnetic screening mass tum was introduced in section© as the magnetic 
analogue of the Debye mass. It gets its leading order contribution from the 
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magnetostatic sector, like the spatial Wilson loop. 

Four dimensional data for SU(2) have been taken j47j . But the numerics 
is much more involved than that for the Wilson loop. Qualitatively the data 
for the screening mass are compatible with the behaviour of the spatial Wilson 
loop tension. In the cold phase its value is about twice the lowest glueball mass. 
Beyond Tc it starts to rise, as you can see from the lower part in fig.©: 



rriM 



rg\T)T 



(7.6) 



as one would expect from a mass in the 3d theory. 

Because the mass is high, the signal to noise ratio becomes small and nu- 
merical extraction becomes tedious. 

The 4D data are being improved [IS] for the region around 2Tc, where we 
want to confront them with the 3d data. 

Quantum numbers of the magnetic screening mass in SU(2) are given by 

= 0^ Precise 3D data for SU(2) are available from ref. |^ and |4Hj . 
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Figure 9: Screening mass as a function of temperature, both in units of Tc, as extracted 
from spatial (top) or temporal (bottom) 't Hooft loops. Below Tc both coincide. The 
arrow gives the mass of the scalar glueball at T = 0. From ref. 

The 4D data are certainly not compatible with the lowest 0^ state. There are 



^^Remember in SU(2) gauge theory all of the gauge-invariant sector has positive charge conjuga- 
tion. This is due to the pseudo reality of SU(2): any element u is equivalent through CT2W(72 = u* 
to its complex conjugate. And the Pauli matrix 1J2 is precisely charge conjugation: A'jf — <T2A^a2- 
As the gauge invariant sector involves integrating over charge conjugation, so is charge conjugation 
invariant. We drop the label for C altogether. 
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c) d) 



Figure 10: (a) is the one loop contribution, and stands for gluon and ghost loop.(b) is 
the two loop contribution. (c) is the renormalization of the thermal Wilson line(double 
circle) inserted (dotted line) into the one loop.(d) is the kinetic term with the loop 
being a gluon or ghost loop 



two recurrences, the highest of which (ref. table 28) is compatible with 
the 4D data at T = 2Tc in fig.®, up to 30%. The situation is currently under 
investigation [S^ . 

7.2 The Debye mass 

The Debye mass to next to leading order is given by: 

mD = {^ + ^)hT + 4^T( log(i±^ + 7.0) + 0{g^) (7.7) 
So 47r g 

The first term is the result from the one loop self-energy discussed in section 
((21). Remember only hard modes contributed to one loop order. 

The second term is only 0{g) smaller. The log term in the coefficient is 
due to Rebhan and comes from scales between electric and magnetic ones. 
The large number comes from all scales equal or smaller than the magnetic 
one, g1^. It is non-perturbative and calculated by numerical simulation on the 
lattice [12] [21] |22] for N=2 and 3. How is explained below. 

A natural choice in 3d is the reduced version of the Wilson line: P = 
TrexpiAf). Its imaginary part is the negative R parity channel. 

• Exercise: Use the electrostatic action jC^; in the correlator. 

Show < ImP{Q)ImP{z) >e has as dominant decay mode exp— 3m£)Z 

The mass can be extracted from this correlator. In terms of the electrostatic 
parameters ea. H6.7() it has a simple form: 

mD = mE{l + dJ- + 0{x^)) (7.8) 
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Figure 11: 't Hooft loop tension pi for SU(2), in units of T^^, as a function of the 
reduced temperature t. The straight hne is a power law fit to t < 1. The fitted 
exponent is 1.32(6), to be compared with 2v ^ 1.26 for the 3D Ising model. The 
curves show the perturbative result, to leading (upper) and next (lower) order, eq. 
()7.21|) in the text. From ref.j 



We have to evaluate the correlator along the physics line eq. (|H.9j) . To see 
that we can expect this form for the correction take the self energy correction 
to one of the three Aq propagators contributing to < ImP{{))ImP{z) >e- Let 
the momentum of this propagator be p. Schematically, the self energy insertion 
reads: 

^^i (27r)3(P+m|)((^+p)2 ^'■'') 

We are interested in the contribution of this integral to the Debye mass in 
the exponent exp— m^z. So we have to evaluate the integral on-mass shell: 
+ mj-) = 0. That contribution is for dimensional reasons ~ That, 

2 J — 

together with the coupling constant factor g"^ gives indeed ^ ~ y'y' 

Here the value of the coefficient d for N=3 is given. Its presence renders the 
next to leading term dominant till T ~ lO^Tc! 

But what is even more so is that the remnant of the series in ea. (|7.7() con- 
verges well up to T = 2Tc. This follows from comparing to the full 4D lattice 
data 137]. 

And so the pattern here is the same as for the spatial Wilson loop: once 
the magnetic sector has contributed to the series the rest is well convergent. 



7.3 Spatial 't Hooft loop 

In this subsection we finally make good our promise to compute the first three 
orders in the tension of the spatial 't Hooft loop at high temperature. We met 
this tension in section (jlj, where it came up as a ratio of twisted partition 
functions, ea. H4.19|l . 
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First we shall give a definition that makes clear its connection with the 
thermal Wilson line and with thermal Z(N) symmetry. 

Consider a closed loop L in the (x,z) plane, and define the spatial 't Hooft 
loop Vfc(L) as a gauge transformation that has a discontinuity Zk on the minimal 
surface spanned by the loop. 

This definition tells us that the loop is a closed magnetic flux loop with 
strength 2;^, as we defined in section (|4.1() . 

A simple realization of such a transformation is the solid angle with 
which L is seen from a point x. It jumps over 47r when crossing the surface. 
So the corresponding operator in Hilbert space is formed by taking Gauss's 
operator: 

Vfc(L) = expi / dx-TrE{x).Du:L{x)^^. (7.10) 
J g 2N 

Remember from section (jlj) that Yfc = diag {k, k, ...,k,k — N, ...k — N) with N-k 
entries k and k entries k-N, so that it generates the center group element: 

27r 27r 
exp (ij^Yk) = ex-vik— = z^. (7.11) 

What does a thermal Wilson line feel when it passes through the minimal 
area of the loop? 

The answer is simple: since the Wilson line represents a heavy test quark 
it will pick up the Z(N) phase Zk- 

So the operator acts like a twist. To recover the twisted partition functions 
flj: of section ^ (with k = (0, k, 0) for the case at hand), just extend S{L) over 
the full X, z cross section As promised in that section we now compute the 
area law for the 't Hooft loop, when T >Tc: 

< Vk{L) >= exp -pk{T)S{L) (7.12) 

and it is the tension Pk{T) we are after. 

7.3.1 The strategy for computing pfc(T) 

Imagine the loop immersed in the plasma. Far away from the loop the value 
of the Wilson line is some fixed Z(N) value, as we learned in section (jSIHl)- We 
will take the value P{Aq) = 1, but any other would have doe equally well. 

The geometry of the problem is such that all profiles orthogonal to the plane 
of the loop will be identical. Only near the border this is no longer true. We 
are only interested in the surface effects, not in the border effects. As we said 
before the Wilson line will jump at the surface of the loop. Such a typical 
profile is shown in fig. (|7.3.1|) . 

The jump should be such, that the derivatives on both sides are equal. Once 
we have fixed these boundary conditions, all we have to do is to determine what 
profile minimizes the effective potential. This is what we are going to do in the 
next section. 



^•^For the vortex correlation in the 2+1 dimensional theory in seetion precisely the same holds 
true. Mutatis mutandis the reasoning of this section applies to that correlation as well. 
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Figure 12: The Wilson line profile, that dominates the steepest descent calculation 
of the effective potential. The parameter q is explained in the text. 

7.3.2 Effective potential for constant profile 

To be specific take SU(2) gauge theory. The key to computing the surface 
tension is the distribution function E{P) of the thermal Wilson line that we 
introduced in section (|3.3j) . eq. (|3.11|) . We are interested in the logarithm of this 
distribution function: 

E{P) = J DA6{Tr exp iC - TrV{Ao)) exp -S{A) = exp -L^{V{q) + /) 

(7.13) 

We have introduced the diagonal traceless matrix C = -7rdiag(g, —q) to 
parametrize P. The free energy density / normalizes the distribution. 
is the spatial volume. The distribution function has two peaks, one at P=l, 
one at P=-l. So the effective action V{q) has minima at q = integer. Even 
(odd) values correspond to P = 1(— 1). 

The constraint in ea. H7.13|) tells us that the potential fluctuates around 
the background matrix CT. So this matrix will be present in the Feynman 
rules as a background field. This is familiar. What is less familiar is that 
the background gets renormalized through the renormalization of the thermal 
Wilson line. This is an effect coming in from two loop order on, and is shown 
in fig. (llUb h Only the sum of the two loop diagrams b) and c) is gauge 
independent. 

The covariant derivative is Dq{Aq = TC). Such a derivative is diagonal in 
the Cartan basis for the fields A = ^ + Q+t' + Q't^ + Q^t^. The Pauh 
matrices r relate to \/2t'^ as ti it iT2. 
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So the Feynman rules consist of the Euchdean ones, except that po is re- 
placed by po ± 2TTTq when it is the Matsubara frequency of Q^. Hence it 
appears in propagators of and in vertices when it pertains to an outgoing 

hue. 

The one loop computation (see fig- HlUb ) is the calculation of the determinant 
of the quadratic part in the quantum fields Q of the action. In Feynman gauge 
you get: 

V{q) + / = log det{-D{Cf). (7.14) 

The choice of gauge does not matter as is shown explicitely in Pisarski's lectures. 
D{C)'^ is the covariant Laplacian. 

There is a well-known identity that relates log det M = Tr log M for any 
matrix M. Then, in the momentum basis you obtain in terms of the Fourier 
transformed Laplacian for V{q) + f: 

/ (^(log((27rr(n-g))2+P)+log((2^r(n+g))2+F)+log((27rrn)2+/^) 

(7.15) 

with the contributions from ,Q integrations written out explicitely. If 
we sum over all values of n, clearly the result will be periodic in q mod 1 and 
even in q. This is of course a consequence of the thermal Z(N) symmetry. In 
fact one gets: 

V{q) = ^7r^TV{l-\q\? (7.16) 

As long as qT = 0{T) the integration over / is hard. For q = 0{g) or g = 
\ — 0{g) we have for n = Oorn = — 1 soft momentum contributions. 

7.3.3 Varying profile and gradient expansion 

Till now we computed the potential as if the profile q were constant. But what 
we are really after is the profile of g as a function of y in between the values 
P = 1 and —1 at y = ±00. That means we have to include the kinetic term 
K = {^dyCy to our potential V to get the tension: 

/oo 
dy{K + V) (7.17) 
-00 

So we tunnel through the potential mountain V{q) from P = 1 to P = —1. 

We swept a little problem under the rug. We forgot the contribution from 
the one loop potential due to the gradient in q{y)\ This however is of higher 
order as we will see shortly. 

To find the tension you have to minimize U with the boundary condition 
that the wall is between P=l and P=-l regions. This is done by the method 
of completing the square: 

U{q) = J dyi{K^/^ - V^/^ f + 2{KV)^/^ (7.18) 
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The integration of the second term over y can be replaced by an integration 
over q, using the chain rule for dyK^/"^ ~ dydyq = Jq dq: 

f dy2{KVfl^ = -(8^22^2)1/2 ^^^1/2 (^^^^g) 

J g JQ 

which gives a number independent of the profile q\ 

Hence U{q) in eq. (|7.18j) is minimized for that profile q that obeys the equa- 
tions of motion 

K^l-i _ = 0. (7.20) 

The tension of the loop is then given by eq. (|7.19|) and using the result (|7.1()j) 
for the potential one gets: 

47r2 

pi(r) = -j=^T\l - C2 as{T) + ...) (7.21) 

Two loop corrections have been computed |38_, from the graphs in fig. HlUb . c 
and d), with the result C2 = 2.0682.... 

What is the typical width of the Wilson line profile? Just look at the 
equations of motion eq. 1)7. 20(1 . Write them out with eq. ()7.16() : 

dyq = inDq{l-\q\) (7.22) 

with m^) = |g'2j'2, the lowest order Debye mass from N = 2. So it is the Debye 
mass that governs the width of the wall. This also answers the question about 
gradient terms. They are there, but are 0{g^) compared to V{q). 

Only hard momenta do contribute to this order, as you can see by computing 
the contribution to pi from the wings of the profile near P = 0, 1. 

The lattice simulated tension |17] in SU(2) is shown in fig. ()ll() . For T ~ 2Tc 
lattice data and perturbative prediction do agree reasonably well. 

The cubic corrections are now known and add a positive contribution. 
They are soft modes on the scale of the Debye mass. As we will see the same 
pattern shows in the pressure, calculated with the same graphs as in fig. HlOf) . but 
with the background q set to zero. The three loop contribution is in progress. 

Magnetic modes will contribute through the next to leading order of the 
Debye mass (see eq. (|7.7|) ) and have not yet been computed. 

For general N and strength k of the loop one finds to one and two loop order 
a remarkably simple scaling law in k |17j : 

Pk{T) = ^ff^Pim (7.23) 

with pi{T) = ^-^^(iV-l)(l- 1.0341. .a^ N + 0{g'^)), periodic in the strength 
k mod N. 

In one loop order this is a simple consequence of additivity of the potentials 
for the various colour modes and of the tunneling path being along the 
direction. Let us look into that in more detail. 
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One tunnels from 1 to by the path qY^, from q=0 to q=l, so the diagonal 
background matrix is C = ^^Yf^. The operator whose determinant we have to 
compute is -D{Cf 

Let us adopt the Cartan basis for the fluctuating potentials Q in NxN matrix 
notation, like we did before for SU(2). This is a colour basis in which the effect 
of the profile is diagonalized in the Laplacian —D'^{C). 

Remember the profile C is diagonal by construction, with diagonal elements 
Cj, i = 1,....,N. The profile appears only in the covariant derivative Dq{C). 
There are diagonal fiuctuations Q'^, that to one loop order do not contribute 
any C dependence: 

Dq{C)Q'^ = doQ'^ + iT[C, Q'^] = dQ'^. (7.24) 

Then there are off-diagonal fluctuations Q*-' (1 <i^j < A^) with Q*-' only 
non-zero in the (ij) entry. For those you have: 

Dii{C)Q'^ = doQ'^ + iT[C, Q'^] = doQ'^ + i{Ci - Cj)Q'^ (7.25) 

and we have diagonalized the Laplacian. The background field C comes in 
through the diagonal elements of the adjoint representation of C = j^Yk. They 
equal 0, or itg up to a factor 27r. 

First we write the contribution to the C dependent part in the potential for 
a fixed combination (ij): 

Vi,{a, -Cj)=TY,J log(r(27rn + Q- Cj)f (7.26) 

The reader will recognize the form of the potential for SU(2), ea. (|7.15|) . 

By changing the sign of the summation over Matsubara modes, it follows 
that Vij{Ci — Cj) is even in Cj — Cj. 

Remember that on the path C = "^jfY^ the profile in its adjoint representa- 
tion equals \Ci — Cj\ = 2Trq , or 0. So the C-dependent part of the full potential 
is obtained by multiplying Vij with the number of off-diagonal modes (ij) with 
a non-zero eigenvalue Q — Cj . 

This number is determined from Yk = diag(fc, ...,k,k — N, ..k — N) by taking 
all differences of elements. There are k elements with value k — N, and N — k 
elements with value k. Hence there are 2k{N — k) ways of picking a non-zero 
combination. We obtain for the full potential: 

V{q) = 2k{N - k){V,j{27rq) - VijiO)) = k{N - k)—T\^{l - \q\f . (7.27) 

The second equality follows from ea. (|7.15|) and ea. (|7.16|) . 

For the kinetic term one has the identity Tr{dyCY = 2^ dy{Ci — Cj)^ , 
valid because TrC = 0. So the same counting as above applies to the kinetic 
term: 

K=- Tr{dyCf = k{N - k)^{dyqf. (7.28) 



*We drop all indices exept colour indices. 
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Then ea . (17.231) follows from the minimization of K + V. 

• The stability group of Yfc is SU{N - k) x SU{k) x U{1). Determine the 
dimensionality of the coset space SU{N)/SU{N - k) x SU{k) x U{1). 
Explain the masses appearing the diagonal and off-diagonal propagators 
-D-\C)). 

Some comments on the 't Hooft loop tension, eq. (|7.23|) 

• The result for the tension is 0{N) in the large N limit, in contrast to the 
tension of the Wilson loop. The latter is 0(1) as follows from the well- 
known index loop counting, valid to all orders in perturbation theory. 

• The absence of a sacred cow is noteworthy in eq. (|7.23|) . For k of order 
1 the result has corrections 0{1/N)\ A priori we would have expected 
corrections of 0{1/N'^) in a theory with only gluons, like in the pressure. 
This anomalous correction is a simple consequence of the counting: not 
N"^ — 1 gluons, but only 2k{N — k) contribute to the tension, whereas they 
all do contribute in equal amount to the pressure. 

Physically ea. (|7.23() is understood as being due to the flux of the screened 
quasi-particles and will be discussed in section Q. 

7.4 Pressure 

In fig. (23} you see the pressure as measured by lattice simulation for three 
colours. It is compared to the analytically computed [20] terms in perturbation 
theory up and including g^log{^). The series starts with the contribution 
from the gluons as free quasi-particles: the Stefan-Boltzmann gas with pressure 
Po = s-s computed from fig. (|10b ) . The interactions between the gluons 

give the two-loop contribution as in fig. HlUb ). The contribution of the Debye- 
screened gluons is the dominant one from /2b in ea. H6.2|) . It equals the logarithm 
of the determinant of the Higgs modes with mass tue- Using dimensional 
regularization 

This being positive and large with respect to the two loop contribution (like 
in the case of the 't Hooft loop) spoils the convergence of the series except for 
academically high temperature. The two next orders g"^ and undo in part 
the effect of the cubic term. But it is clear that perturbation theory has little 
predictive power in the region 2Tc to lOOTc, as it varies there over more than 
20%! 

To put the calculation of the contributions of order higher than three in 
perspective and to see how the different scales come in, we recall once more 
the hierarchy of scales, cut-offs A and reduced actions needed to compute the 
pressure: 

T » Ae » gT » Am >> g^T 
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Figure 13: Left: perturbative results at various orders, including 0{g^) for an optimal 
constant. Right: the dependence of the 0{g^) result on the (not yet computed) 
constant, which contains both perturbative and non-perturbative contributions. The 
4d lattice results are from . From ref . [201 ■ 



The pressure is normahzed by po and consists of three parts: 

P 

— = Ph+PE +PM 

Po 

The hard modes are cut-off in the infrared by A^; and equal ph- SchematicaUy 
we get: 

T T 
Ph = l + g'^ + g'^log — + g^ + gHog— + g^ + .. 

All powers of the coupling are even. The short distance scales (larger than 
T) are absorbed in the running coupling, eq. (|6.4jl . The cut-off A^; appears only 
in logarithms. The electric mode contributions are computed with Ce and give 

PE- 

Pi. = / + /log^+5^+/ + /log^ + /log^ + / + .. 
niE rriE Am 

The dominant cubic term was computed in eq. 1)7. 29(1 . We can expect logarithms 
of the two ratios of the three scales itle, ^e and Am in the electrostatic action. 
Finally the magnetic contribution is computed with Cm- 

PM = g log — -Kg" + ... 
9m 
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We only put in the obvious dependence on the parameters in the electro- 
static and magnetostatic actions. There are three comments: 

• All terms shown are perturbatively calculable, except the last one in pM- 

• All perturbatively calculable terms have been computed [301, except for 
the terms. In particular the log's are known by now |2()j . 

• All dependence on the cut-offs cancels, as expected. 
Clearly this is quite a calculational performance! 

So at this point one expects that the miracle of the Debye mass may ma- 
terialize: compute the non-perturbative term of 0{g^) with Cm on the lattice, 
and the perturbative term of 0{g^). Does the ensuing series stabilize down to 

To see whether this may work at all, values for the sum of the yet unknown 
non-perturbative and perturbative coefficients have been put in the right hand 
figure H13|). Clearly there is a window where lattice data do connect smoothly 
to the series. 



8 Flux of quasi-particles as seen by spatial 
Wilson and 't Hooft loops 

The idea of quasi-particles is to leading order in the pressure embodied by the 
Stefan-Boltzmann term. It counts all the degrees of freedom of the gluons, 
colour and spin. 

What we want to argue below is that specific degrees of freedom of the 
gluons, namely their flux, are clearly seen in the behaviour of the spatial 't 
Hooft loop, especially the scaling law in the strength of the loop. Then we will 
look for a similar effect in the Wilson loop, but now with quasi particles that 
carry magnetic flux. 



8.1 Gluon flux and the 't Hooft loop scaling law 

We found in perturbation theory including two loop order: 

P.(r) = ^^f^Pi(r). (8.1) 

Below we will show that the factor A;(A^ — A;) is the number of gluons in the 
adjoint multiplet with a fixed value N for the charge characterizing the strength 
k of the loop. Essential is the physical meaning of the loop: it measures the 
colour electric fiux in the plasma. The loop, ea. (|7.1U(l . can be rewritten as a 
sheet of electric dipoles: 

VkiO = ^ / dS.TrEYk (8.2) 
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and measures the colour electric flux going through the loop. 

Expression ()8.2() is a dual Stokes law: remember the original definition of the 
't Hooft loop is that of a loop of a Dirac magnetic flux. But there is no magic 
about (|8.2jl ! Both definitions can be shown to have the same commutation 
relation with the Wilson loop 

Vk{C)W{C')Vk{C)^ = zl^^'^"^W(C') (8.3) 

n(C, C) being the number of times C and C loop each other. 

Hence the product Vk(C)V\C^ does commute with the Wilson loop. 

So the product is a regular gauge transformation and acts as the unit op- 
erator in the physical Hilbert space. Both the representation ea. ()7.1U|) and the 
flux representation ea. (|8.2() are identical there. 

• Show that Vfc(C) has the 't Hooft commutation relation ea. (|8.3() with a 
spatial Wilson loop W(C'\ (Hint: use the canonical commuation relations 
between field strength E and potential A in the Wilson loop.) 

Consider a gluon close to the minimal area of the loop. As its flux is screened 
it needs to be within the screening length Id to shine its flux through the loop 
The gluon is in the adjoint representation so its charge 1^ is either or itA'^ as 
you can see from the differences of eigenvalues of . There are then obviously 
2k{N — k) gluon species in the adjoint representation with this charge itA^. 
Each one individually shines flux through the loop. This flux equals zizN/2, the 
other half of the flux is lost on the loop. That means that such a gluon will 
give a contribution to the loop of 

Vk(.C)\one gluon = exp((i^(±y)) = -1. (8.4) 

Assume that the distribution P(l) of gluons in the slab of thickness Id 
containing the loop is independent of the species and say Poissonian: 

P(/) = ir'exp-[. 

Here / is the average number of gluons of that species in the slab. Then the 
average of the loop over just one species is: 

< Vfc(C) > lone species = ^J^l'^exp-l = exp -2/". (8.5) 

If the contribution of all 2k{N—k) species is independent the result becomes: 

< Vk{C) >= exp -4k{N - k)l. (8.6) 

^^This is a simplification. Also gluons farther out can contribute, giving a change in the overall 
factor of our final result. 
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Now / = 2l£)S{C)n[T). n(T) is the density of a specific gluon species and an 
area law results with the tension: 

Pk{T) = 8k{N - k)lDn{T) (8.7) 

This formula represents the physical raison d'etre of the tension. It is due to 
the screening of the electric flux of the gluons, their density and their degeneracy 
with respect to the charge characterizing the strength of the loop. The 
Poisson distribution function is not essential to this result. Any thermodynamic 
distribution function, that is peaked around the average T with a width — (J)"^ 
proportional to will give an area law. It is this proportionality constant that 
will appear in ea. (|8.7() . But it does not depend on the strength k and will drop 
out in ratios. 

Hard gluons with momentum of order T, will have a photon like distribution 
function 

p(/)~(/7(i+oy. 

As the reader can easily check, it has a variance P — = so fluctuations 

of order 1. Not surprisingly it does not give an area law as you can find out by 
plugging it in ea. (|8.5() . Small fluctuations are essential for the area law. 



8.2 Magnetic flux and k-scaling of the spatial Wilson 
loop 

The question is now: what about similar ratios of the spatial Wilson loop? 

We convinced ourselves in section (fTj) that the leading contribution to its 
tension came from the static magnetic sector. This sector is populated with 
static transverse gluons, i.e. gluons with a static magnetic field screened in 
some non-perturbative way with a screening mass mjv/- The screening mass 
is inversely proportional to the coupling gf .j on dimensional grounds. More 
precisely, in our discussion of the magnetic screening in section ((T)) we noticed 
that the proportionality constant was large with respect to other scalar masses, 
because of its unnatural parity. 

In the context of our model for the 't Hooft loop tension, it is interesting 
to see what the ratio of the tension of the Wilson loop to magnetic screening 
mass is according to the lattice 

^ = 8.15(15) SU(2) (8.8) 

Suppose that the screening length Im = rn'^j of our magnetic quasi-particles 
is much smaller than the mean distance between them: 

nulli <- 1- (8-9) 

Both are parametrically equal to g\j. There must be a dynamical reason for 
this ratio to be small. It supposes the magnetic screening of the quasi-particles 
is so efficient that they constitute a gas of approximately free lumps. 
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If the latter is true, then the reasoning in the previous section would, mutatis 
mutandis, give for the tension of the Wilon loop: 

a = c hiUMiT) (8.10) 

where c is some numerical constant depending on your preferred probability 
distribution. So this relation turns the numerical result for the screening length 
eq. (|8.8j) into a numerical result for the density of screened magnetic quasi- 
particles: 

c tom/m = 0.015.... (8.11) 

This means one quasi-particle in about seventy screening volumes if c=l, an a 
posteriori justification for the model. For the Poison distribution this density 
is even smaller. 

This is for SU(2). 

Remember that for SU(3) and higher groups the magnetic screening was 
related to states. From Teper's work |46| the lowest 0^^ mass in units of 
the string tension for SU(N) [N > 3)) gives: 

umIm = 0.028.. ..(1 - (1.64...)/iV2 ^ o{N-^)) (8.12) 

Although both screening mass and square root of string tension are para- 
metrically the same, their ratio is for dynamical reasons small! So we may 
say that the quasi-particles are to a good approximation free. How good will 
become more clear at the end of this section. 

So having taken courage we now turn to a flux representation of the spatial 
Wilson loop jl2] In fact, if we want to find out about the strength or N- 
allity k of the loop, we would expect in analogy with the 't Hooft loop, ea. ()8.2|) . 
a magnetic dipole distribution projected on the k-charge: 

Wk{C)^expi^[ dS.TrBYk (8.13) 

iV Js{C) 

with the magnetic field strength B replacing its electric counterpart E. This 
is incorrect for at least one reason: as discussed in section the Wilsonloop 
with N-allity k defines a periodic tension only after averaging. Now the left 
hand side of the equation is periodic. So the equation must be considered in 
thermal-averaged form. There is much more to say about this formula! For a 
derivation of a slightly different form see ref . • 

But if we accept this flux representation on the basis of its analogy with the 
't Hooft loop, then the same reasoning as for the 't Hooft loop applies and one 
concludes that the tension ak scales as its electric analogue: 

<^k{T) = J^^'^iiT) (8.14) 

All we have to do now is to hear the verdict of the lattice; the ratios found 
by simulation |^ are close- within a percent for the central value-: 

5[/(4) : £72/^1 = 1.3548 ± 0.0064 
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SU{6) : = 1.6160 ± 0.0086; a^/ai = 1.808 ± 0.025 

The results are that precise, that you see a two standard deviation, except 
for the second ratio of SU(6). As we said, magnetic quasi-particles are dilute 
but only approximately free. 

There is a less precise determination of the ratio o"2/ci = 1-52 it 0.15 in 
SU(5) jSB]- But the central value is within 1 to 2% of the predicted value 3/2. 
In conclusion, lattice results are quite encouraging. More simulations are under 
way [52| . 

9 Epilogue 

We highlighted some of the salient problems in the critical region of deconfined 
and chirally restored QCD. Universality works well, when it is defined at all. 

In the region above 2Tc we studied the indications that perturbation theory 
works well if the excitations of the magnetic sector are included. These excita- 
tions, the magnetic quasi-particles ("magnons"), may be considered as screened 
quasi-particles. 3d simulations suggest that their screening is so strong that 
they form approximately independent lumps. This idea is strikingly vindicated 
by the dynamically small ratio of magnetic screening over string tension. 

The idea works quite well in the Wilson loop sector. And the result for the 
Debye mass again supports this idea; including the magnetic quasi-particles cor- 
rects a result which was off by more than 0(1) without! The non-perturbative 
part of the pressure to order is the partial pressure of the magnetic gluon 
gas. 

Unfortunately our arguments are at present only good enough for parametric 
statements. In cases like the Wilson loops we get rid of the unknown constants 
through ratios. A curious situation prevails: we are indirectly witness of the 
magnetic quasi-particles, but they remain for the time being the Poltergeists of 
the magnetic sector. Of course they are only approximately free: their density 
and pressure are subject to a logarithmic rescaling of the temperature g'^{T)T, 
and hence at asymptotic temperatures the Stefan Boltzmann limit prevails. 

10 Suggestions for further reading 

Textbooks: 

• J. Kapusta, Finite Temperature Field Theory, CUP 1989. 

• M. Lebellac, Thermal field Theory, CUP 1996. 
Schools: 
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• M. Shaposhnikov (Erice 1996), High Temperature Effective Field Theory. 

• M. Laine (Trieste 2002), Finite Temperature Field Theory. 

• F. Karsch, Lect. Notes Phys. 583 (2002) 209-249, Lattice QCD at High 
Temperature and Density. 



K.Kanaya, An introduction to finite temperature QCD on the lattice. 



Prog.Theor.Phys.Suppl.131; 73, 1998; |hep-lat7 9804006| 



• Robert D. Pisarski, Notes on the Deconfining Phasetransition. Proceed- 
ings of Cargese Summer School on QCD Perspectives on Hot and Dense 
Matter, Cargese, France, 6-18 Aug 2001; hep-ph/0203271 

• A.K. Rebhan, Hard thermal loops and QCD thermodynamics. Proceed- 
ings of Cargese Summer School on QCD Perspectives on Hot and Dense 
Matter, Cargese, France, 6-18 Aug 2001; hep-ph/0 111341| 
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Appendix A: Free energy of heavy source 
and Wilson line average 

Below we prove eq. (|3.9|) in the main text. It relates the free energy of the gauge 
averaged heavy quark state to the thermal average of the Wilson line. 

Start with the gauge average of a state with one heavy quark. The state 
with one quark, and eigenstate of the field operator A is 

4\A). (10.1) 

We have suppressed reference to the position x of the quark. Gauss' opera- 
tor is G(^o) = TrE.DAo + iI^^AqiP for an infinitesimal transformation. Aq is 
dimensionless and is up to a dimensionful multiplicative constant the scalar 
potential at Euclidean time r, as in eq. (|3.7j) . Integration over three-space is 
understood. To render the state, ea. HlO.lj) . gauge invariant, we have to take 



53 



the gauge projector Pq = J DAQexpiG{AQ) and act with it on the state. So 
the free energy excess due to the quark becomes: 

1 ^ 

exp -F^ = -J2 TT^{A\^iPG exp -H/TPg4\^) /Tt^{A\Pg exp -H/TPg\A) 

(10.2) 

The gauge projector and the Hamiltonian commute. As in the usual tran- 
scription of the free energy into path integral language, we write for a partition 
of the Euclidean time interval into n bits <5t = 

Pg exp —H/TPg = Pg exp —StHPg exp —6tHPg Pg exp —StHPg 

(10.3) 

using that Pg is a projector. The transcription is identical to that of the free 
energy except that after rewriting the n Hamiltonian factors as a path integral 
we are still left with the n fermionic projectors ex.p iip"^ AqiP . Every single one 
of them has the effect - through the canonical anticommutation relations-: 

exp {ixl:Uoi^)4\A) = i;l{expiAo)ki\A). (10.4) 

The reader will recognize the string bit appearing in ea. H3.7|) Repeat the op- 
eration in eq.(^I31) n times and the net result is the familiar time ordered 
exponential ^^{Ao) multiplied on the left by the fermion operators: 

A^i{ViAo))ki. 

Using once more the anti-commutation relations, one finds the anticipated rela- 
tion between the thermal average of the Wilson line and the free energy excess: 

exp-F^ = J DAexp-^SiA)^TTV{Ao)/ J DAeicp-^SiA) = {P{Ao)). 

(10.5) 

as in the main text, eq. (|3.9j) . 
Corollary 

What is true for one heavy quark can be shown along the same lines to be true 
for the correlator of heavy quarks: it is given by the correlator of two Wilson 
lines. 
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